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Preface 


The  purpose  of  this  study  was  to  develop  and  test  an 
optical  position,  scale,  and  rotation  invariant  pattern 
recognition  system  using  a  computer  generated  hologram  and 
spatial  light  modulator.  This  research  contributes  to  the 
Air  Force's  on-going  investigation  of  autonomous  target 
recognition  systems. 

I  received  a  tremendous  amount  of  help  and  encourage¬ 
ment  while  working  on  this  thesis.  First,  I  am  deeply 
indebted  to  my  thesis  advisor.  Dr.  Matthew  Kabrisky,  for  his 
patience  and  understanding  help  with  all  aspects  of  the 
research.  I  would  also  like  to  thank  the  members  of  my 
committee.  Dr.  Steven  K.  Rogers  and  Dr.  James  P.  Mills,  for 
their  guidance  and  timely  suggestions,  especially  with  the 
optics.  My  gratitude  goes  to  Captain  Robert  Williams  for 
his  insightful  comments  and  extended  help  with  the  digital 
design  and  modifications  of  the  Litton  spatial  light 
modulator.  I  would  also  like  to  thank  Captain  Dave  Miazza 
for  his  guidance  and  help  with  the  overall  project. 

Finally,  and  most  importantly,  I  would  like  to  thank  my 
wife,  Carey,  for  her  understanding,  encouragement,  and  help 
during  this  thesis  effort. 


li 

•Original  conl.uiuG  color 
platos:  All  DTIC  reproduct¬ 
ions  will  bo  in  black  and 
white" 


Table  of  Contents 


page 


Preface  .  ii 

List  of  Figures  .  v 

List  of  Tables  .  vii 

Abstract  .  viii 

I.  Introduction  .  1 

Sununary  of  Current  Knowledge  .  2 

Problem  Definition  .  5 

Scope  of  Thesis  .  5 

Approach  .  6 

II.  Operation  of  the  SLM .  8 

Construction  and  Basic  Operation 

Principles  .  8 

Modification  of  the  SLM  .  10 

Litton  Signals  .  13 

Interface  with  the  Semetex  board  .  14 

III.  Polar-Logarithmic  Coordinate  Transformation  CGH 

for  Scale  and  Rotation  Invariance .  19 

Design  of  the  lnr-9  Coordinate 

Transformation  CGH  .  19 

Experiment  and  Results  .  21 

Conclusion  .  30 

IV.  Position,  Scale,  and  Rotation  Invariant 

Feature  Space  .  31 

Introduction  .  31 

Position  Invariance  .  33 

Scale  and  Rotation  Invariance  .  33 

Optical  Processor  for  PSRI 

Feature  Space  .  37 

Experimenatal  Results  .  39 

Conclusion  .  48 

V.  Optical  Correlation  .  49 

Introduction  .  49 

Theory  .  50 

Experimental  Configuration  for  Optical 
Correlator 


53 


VI.  Conclusions  and  Recommandations  . 

Conclusions  . 

Recommendations  . . 

Appendix  A:  Software  Development  for  Litton  SLM 

Appendix  B:  Fabrication  of  the  CGH  . 

Bibliography  . 

Vita  . 


l' 


> 

»fj 

i'i 

♦MJ 


Figure 


List  of  Figures 


VanderLugt  correlator . 

Operation  of  LIGHT  MOD  as  SLM . 

Block  diagram  of  128x128  LIGHT  MOD . 

LIGHT  MOD  timing  diagram . 

Semetex  interface  board  schematic  diagram. . 

Pixel  write/erase  control  through  Semetex 
interface  board  to  LIGHT  MOD  device . 

Lnr-d  CT  realization . 

Scale  and  rotation  invariant  setup . 

(a)  Input  square  (b)  167%  scale  change 
and  45*  rotation . 

Theoretical  Inr-G  CT  of  squares . 

Experimental  lnr-0  CT  of  squares . 

(a)  Input  letter  X  (b)  45*  rotation  of  X... 

Theoretical  Inr-G  CT  of  X’s . 

Experimental  lnr-9  CT  of  X’s . 

PSRI  feature  space  optical  processor . 

(a)  Input  square  (b)  Scaled  and  rotated 
square . 

Fourier  transforms  (a)  square  (b)  scaled 
rotated  square . 

PSRI  feature  space  (a)  sqaure  (b)  scaled 
rotated  square . 

(a)  Template  of  F15  (b)  Scaled  and  rotated 

version  of  template . 

F15  scene . 

PSRI  feature  spaces  (a)  F15  template 
scaled  and  rotated  F15  (c)  real  scene . 

Optical  correlator  using  MSF. . 


5.2  MSF  construction  using  Vander  Lugt  filter 


51 


5.3  PSRI  feature  optical  processor  and  optical 


correlator  using  MSF's .  55 

B.l  Laser  printer  plot  of  the  CGH .  82 

B.2  Optical  reduction  of  CGH .  83 


List  of  Tables 


Table 


2.1  Required  signals  from  computer  interface  to 

pin  J5  on  the  X  Driver  Card .  12 

2.2  Cable  assembly  Cl.  Semetex  interface  card  to 

X  driver  card  .  18 


Abstract 


This  thesis  investigates  the  integration  of  an  optical 
system  for  real-time  position,  scale,  and  rotation  invariant 
pattern  recognition.  Specifically  a  Litton  Magneto-Optic 
Spatial  Light  Modulator  is  interfaced  to  a  Zenith  248 
microcomputer  and  AT&T  frame  grabber.  A  user  interface 
written  in  C  allowed  arbitrary  patterns  to  be  written  to  the 
SLM  by  a  user.  Patterns  can  be  generated  in  computer  code, 
transferred  from  a  CCD  camera  (via  the  frame  grabber) ,  taken 
from  an  image  file  on  computer  disk,  or  taken  from  a  VAX 
image  file.  The  ability  to  locally  generate  a  computer 
generated  hologram  using  th«  computer  generated  interfero- 
gram  method  on  the  Sun  workstations  and  Imagin  laser  printer 
was  developed.  The  CGH  developed  was  used  in  an  optical 
system  which  performed  a  Fourier  transform  and  lnr-& 
coordinate  transformation  to  create  a  position,  scale,  and 
rotation  invariant  feature  space.  CCD  cameras  were  used 
throughout  the  optical  processor  for  display/analysis .  The 
entire  system  was  developed  and  tested  using  various 
templates  and  real  scenes. 
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COMPUTER  GENERATED  HOLOGRAM  AND  SPATIAL  LIGHT  MODULATOR 
FOR  OPTICAL  PATTERN  RECOGNITION 

I .  Introduction 

Automatic,  real-time  pattern  recognition  is  a  major 
concern  of  the  Air  Force  in  the  area  of  pilot  aids,  smart 
weapons,  and  general  purpose  robotics.  A  pattern 
recognition  system  that  works  in  real-time  and  operates 
independently  of  the  position,  scale,  and  in-plane  rotation 
of  the  input  target  is  needed  to  accurately  find  and 
identify  targets  embedded  in  their  natural  surroundings. 
Pattern  recognition  to  date  has  been  implemented  primarily 
on  large  computer  systems  employing  algorithms  which  require 
an  enormous  amount  of  calculations,  leading  to  a  time- 
consuming  task.  Pattern  recognition,  performed  by  a 
coherent  optical  system,  allows  for  parallel  processing  and 
high  speed  calculations  (nanoseconds) .  These  advantages 
over  electronic  serial  processing  techniques  make  the  use 
of  optics  appealing  for  the  pattern  recognition  task. 

This  thesis  documents  an  effort  to  implement  portions 
of  a  position,  scale,  and  rotation  invariant,  real-time, 
optical  pattern  recognition  system  using  computer-generated 
holograms  (CGH)  and  a  magneto-optic  spatial  light  modulator 


(SLM) .  The  SLM  will  be  used  to  input  real-time  scenes  into 
the  system,  while  the  CGH,  along  with  a  Fourier  transform 
(FT)  lens,  will  perform  the  lnr-0  coordinate  transformation. 
The  general  procedure  to  be  followed  in  solving  the  optical 
pattern  recognition  problem  will  be  discussed  in  the 
Approach  section  of  this  thesis. 

Summary  of  Current  Knowledge 

Many  pattern  recognition  systems  (both  digital  elec¬ 
tronic  and  optical)  rely  on  correlation  using  a  template 
matching  scheme  to  "recognize”  a  given  input.  Correlators 
have  the  inherent  ability  to  accomodate  translational 
shifts,  and  to  simultaneously  handle  multiple  targets.  A 
high  intensity  pea)c  is  produced  during  autocorrelations 
(correlation  with  identical  target)  whereas  low  intensity 
distributions  arise  during  cross-correlations  (correlation 
with  different  targets).  However,  correlators  are  very 
sensitive  to  changes  or  distortions  between  the  ideal 
reference  and  the  actual  input  scene. 

The  optical  correlation  of  two  2-dimensional  images  is 
achieved  in  parallel  and  real-time  by  placing  the  conjugate 
Fourier  transforms  of  reference  objects  into  the  transform 
plane  of  an  optical  processor  like  the  one  shown  in  Figure 
1.1.  In  this  system,  the  Fourier  transform  of  the  input 
scene  is  multiplied  by  the  conjugate  FT  of  the  reference 
object.  The  correlation  of  the  two  space  functions  appears 
in  the  output  plane  (2:027).  The  use  of  spatial  light 


modulators  (SLM)  and  liquid  crystal  televisions  (LCTV) 
allows  a  system  to  handle  real-time  input  scenes  while 
performing  pattern  recognition  with  multiple  reference 
objects  or  targets  (8:466,  12:698-704,  11:467-468). 

The  conventional  VanderLugt-type  correlator  used  for 
pattern  recognition  is  unable  to  recognize  scaled  or  rotated 
Images  of  a  reference  object.  It  has  been  shown  that  a  1% 
scale  change  of  the  reference  object  causes  the  signal  to 
noise  ratio  of  the  resultant  correlation  peak  to  decrease  by 
lOdB  from  that  of  the  autocorrelation.  Likewise,  a  20dB 
loss  occurs  for  a  2  rotation  of  the  input  from  the 
reference  (2:1652).  One  solution  is  to  perform  a  Inr-^ 
coordinate  transformation  preprocessing  operation,  which 
results  in  partial  scale  and  rotation  invariance  (5:289). 
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Figure  1.1.  VanderLugt  correlator 

A  low  cost,  commercially  availible  LCTV  has  been  used 
as  a  SLM  in  a  VanderLugt-type  coherent  optical  correlator 
(Figure  1.1).  Results  have  shown  that  the  LCTV  resolution. 


vK* 

vS' 

J’S* 

V/i" 


'’SSa 

m 

M 

w 

*  i 
^12:1 


’S*S 


contrast,  and  speed  are  sufficient  for  basic  (geometrical 
objects  in  low  noise  fields)  real-time  pattern  recognition 
applications  (11:467).  The  LCTV's  main  advantage  is  its 
versatility.  It  can  be  used  in  the  input  plane  and/or  in 
the  transform  plane,  and  it  can  be  addressed  from  a  remote 
video  camera,  a  personal  computer,  or  an  external 
transmitter  via  a  standard  antenna.  Disadvantages  of  using 
the  LCTV  in  high  resolution  processes  include:  (1)  the 
inability  to  obtain  more  than  a  lOdB  dynamic  range  (12:41), 
(2)  the  Introduction  of  phase  distortion  in  the  passing 
wavefronts  (6:938),  and  (3)  its  sensitivity  to  the 
polarization  angle  of  the  incident  light  (12:72). 

The  Litton  iron  garnet,  H-triggered,  magneto-optic 
device  (LIGHT-MOD'")  has  been  used  in  similar  capacities  in 
more  expensive  optical  signal  processing  systems  where 
greater  dynamic  range,  resolution,  and  speed  are  necessary 
(8:466,  12:698-704,  16:55-64).  The  LIGHT  MOD  functions 
predominantly  as  a  spatial  light  modulator  (SLM) .  The  data 
(either  scene  or  filter)  are  transferred  from  a  personal 
computer  to  the  device  through  the  horizontal  and  vertical 
electrode  structure  of  the  LIGHT  MOD.  Since  each  pixel  or 
line  can  be  switched  in  approximately  Ifisec,  a  512x512  frame 
can  be  transferred  to  the  device  in  0.512  msec  in  parallel 
(17:58).  By  using  the  parallel  addressing  scheme,  scenes  can 
be  input  quic)cly  to  the  processor  when  the  LIGHT  MOD  is 


positioned  at  the  input  plane.  Psaltis  used  two  SLM's  in  a 


VanderLugt  correlating  configuration,  which  illustrates  the 
usefulness  of  the  LIGHT  MOD  functioning  as  a  SLM.  The  input 
image  was  recorded  on  a  SLM  placed  in  the  image  plane  of 
Figure  1.1,  whereas  the  FT  of  a  reference  image  (acting  as  a 
matched  filter)  was  recorded  on  the  second  SLM  placed  in  the 
transform  plane  (16:699). 

The  recent  use  of  computer  generated  holograms  (CGH's) 
in  optical  processors  has  allowed  systems  to  perform 
coordinate  transformations  in  nanoseconds  as  opposed  to  the 
required  calculations  taking  tens  and  hundreds  of  seconds  on 
electronic  digital  computers  (3:217-222,  6:938-942,  17:8). 
Computer  generated  phase  hologram  masks  with  a  given  phase 
transfer  function  can  be  used  to  produce  various  geometri¬ 
cal  transformations.  The  desired  coordinate  transformation 
for  this  thesis  was  the  Inr-G. 

Problem  Definition 

The  thesis  problem  is  the  development  of  portions  of  an 
optical  pattern  recognition  system  which  can  locate  a  given 
target  independent  of  its  position,  scale,  or  rotation  in 
realistic  input  scenes.  The  system  must  be  able  to  recog¬ 
nize  and  locate  the  target  in  real-time.  The  real-time  goal 
requires  integration  of  a  CGH  and  a  computer  controlled  SLM. 

Scope  of  Thesis 

This  thesis  focused  on  constructing  a  real-time, 
position,  scale,  and  rotation  invariant  (PSRI)  optical 


pattern  recognition  system  to  identify  arbitrary  targets  in 
low  noise  scenes.  The  integration  of  a  CGH  into  an  optical 
system  and  the  control  of  a  SLM  by  a  computer  was  a  main 
thrust  in  constructing  the  real-time  system.  In  general,  a 
generic  optical  processing  lab  was  established  to  carry-out 
the  ongoing  pattern  recognition  effort  at  AFIT. 

Approach 

The  thesis  problem  was  split  into  three  different 
parts.  The  first  part  included  the  LIGHT  MOD  modifications, 
software  development  for  the  LIGHT  MOD,  and  an  investigation 
of  using  the  LIGHT  MOD  as  an  input  and  transform  plane 
device.  Particular  emphasis  was  placed  on  integrating  the 
LIGHT  MOD  into  the  Zenith  248/AT&T  frame  grabber/Sony  CCD 
camera  system  and  creating  a  user  friendly  environment. 

Part  two  concentrated  on  reconstructing  a  rotation  and 
scale  (but  not  position)  invariant  optical  system  as 
documented  in  the  literature  (4;  5;  13) .  The  development 
and  reduction  of  a  lnr-0  CT  CGH  was  included.  Results  were 
compared  to  those  obtained  by  others  in  the  literature  to 
verify  the  validity  of  the  CGH. 

Part  three  was  the  PSRI  optical  system  implementation. 
Targets  were  placed  into  the  system  as  transparencies  or 
cutouts.  The  magnitude  squared  of  the  Fourier  transform 
of  the  input  image  was  detected  by  a  CCD  camera  and  sent  to 
the  LIGHT  MOD  which  was  functioning  as  a  SLM.  The  magnitude 
of  the  Fourier  transform  was  then  used  as  the  input  to  a 


lnr-&  coordinate  transformation  computer  generated  hologram 
in  order  to  construct  the  position,  scale  and  rotation 
invariant  feature  space.  This  feature  space  could  then  be 
used  as  the  input  to  an  optical  correlator.  An  optical  cor¬ 
relator  was  constructed  during  the  thesis  effort  but  no 
correlating  results  were  obtained  due  to  insufficient  time. 


In  this  chapter,  properties  and  modifications  of  the 
Litton  iron-garnet,  H-triggered  magneto-optic  device 
(LIGHT-MOD)  are  characterized.  In  the  first  section,  the 
construction  and  operation  principles  of  the  LIGHT  MOD  will 
be  discussed.  The  second  section  addresses  the  modifica¬ 
tions  required  to  use  the  LIGHT  MOD  with  the  Zenith  248 
computer  system.  The  last  section  briefly  discusses  the 
performance  of  the  LIGHT  MOD. 

The  LIGHT  MOD  can  be  used  in  optical  signal  processing 
as  an  electrically  alterable  high  speed  two-dimensional 
spatial  light  modulator  (SLM) .  Such  a  SLM  can  function  as 
an  electrical-optical  interface  to  input  Images  into  the 
object  plane  of  an  optical  system,  to  place  spatial  filters 
into  the  Fourier  transform  plane,  and  to  modulate  or  scan 
the  processed  output  image  or  wavefront.  The  LIGHT  MOD  in 
this  thesis  effort  functioned  as  a  SLM  to  input  real-time 
scenes  as  well  as  to  strip  the  phase  from  the  Fourier 
transform  of  input  scenes.  The  magnitude  squared  of  the 
Fourier  transform  was  displayed  on  the  SLM  and  used  as  the 
input  to  the  lnr-9  coordinate  transform  CGH  (CT  CGH) . 

Construction  and  Basic  Operation  Principles 

The  LIGHT  MOD  consists  of  a  bismuth  substituted 
transparent  iron  garnet  film  grown  on  a  nonmagnetic 
substrate.  The  direction  of  uniaxial  anisotropy  is  oriented 
perpendicular  to  the  plane  of  the  film.  Pixels  are  set  in 
even  rows  and  columns  (128  rows  by  128  columns  -  16,384 


total  pixala)  and  metallic  gold  conductors  are  deposited 
along  one  edge  covering  the  full  length  of  the  roes  and 
columns.  The  conductors  cross  at  one  corner  of  each  pixel 
where  a  small  loop  in  the  conductor  increases  the  magnetic 
flux  due  to  current  flow  in  the  conductor  (17:  55-56). 

Operation  of  the  LIGHT  MOD  when  used  as  an  SLM  or  light 
valve  is  depicted  in  Figure  2-1.  A  source  of  polarized 
light  is  required  as  an  input  and  an  analyzer  is  placed  on 
the  opposite  side  of  the  light  source.  As  the  polarized 
light  passes  through  the  pixels,  the  polarization  of  the 
transmitted  light  is  rotated  clockwise  or  counterclockwise 
depending  upon  the  magnetization  within  the  pixel.  The 
amount  of  the  resulting  rotation  depends  upon  the  Faraday 
constant  of  the  material  and  the  thickness  of  the  film 
(17:58).  The  analyzer  is  rotated  with  respect  to  the 
polarization  of  the  input  (laser  source)  so  that  light  will 


Figure  2.1.  Operation  of  LIGHT  MOD  as  SLM  (17:  57) 


either  pass  through  or  be  blocked,  depending  on  the  magnetic 
polarity  of  any  given  pixel.  The  analyzer  is  adjusted 
to  set  image  polarity  (light  scene  on  dark  background  or 
vice  versa)  and  contrast  ratio  (dynamic  range) . 

A  pixel  is  changed  by  simultaneously  applying  a  current 
pulse  to  a  selected  x  and  y  drive  line.  This  current  pulse 
nucleates  or  partial  changes  the  selected  pixel  into  a 
"neutral"  state.  This  state  must  be  followed  by  the  appli¬ 
cation  of  an  external  magnetic  field  (coil  pulse)  to 
complete  the  change  of  the  pixel  to  the  on  or  off  state. 

The  application  of  the  coil  pulse  is  referred  to  as  pixel 
saturation.  The  switched  pixel  will  then  remain  in  that 
state  Indefinitely  until  it  is  nucleated  in  the  opposite 
direction  by  drive  line  currents. 

Modification  of  the  SLM 

Although  the  LIGHT  MOD  was  to  be  a  fairly  small  part  of 
the  optical  pre-processing/pattern  recognition  system,  the 
modifications  required  to  get  the  LIGHT  MOD  operational  and 
user  friendly  became  a  very  large  part  of  this  thesis 
effort.  The  L135,  128x128  display  LIGHT  MOD  purchased  from 
Litton  Data  Systems  is  a  developmental  model  designed  for 
the  Apple  computer;  however,  the  AFIT  labs  are  equipped  with 
Zenith  248  computers  which  are  IBM  AT  compatible.  Available 
to  us  was  the  interface  board  from  Semetex  Corporation.  They 
make  a  similar  SLM  device  which  is  driven  by  the  IBM  class 
computers;  however,  the  Semetex  interface  board  and  the 
LIGHT  MOD  were  not  compatible.  Instead  of  designing  an 


interface  board  to  link  the  Zenith  248  and  LIGHT  MOD,  the 
existing  outputs  from  the  Semetex  board  were  modified  and 
used  as  the  inputs  to  the  LIGHT  MOD  driver  cards. 
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A  block  diagram  of  the  128x128  LIGHT  MOD  display  system 
is  shown  in  Figure  2.2.  Cable  Cl  connects  the  display 
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pulse  (write  pulse,  PT) ,  pixel  saturation  pulse  (coil  pulse, 
PC),  and  grounds.  Table  2.1  shows  the  required  signals  for 
the  X  driver  card  and  their  respective  locations  on  pin  J5 
(a  26  pin  connector) . 


Table  2.1 

Required  signals  from  computer 
interface  to  pin  J5  on  X  Driver  Card 


J5  pin 
number 


Signal 

name 

Cable  signal 
function 

FPBO 

X  Address 

LSB 

FPBl 

•f 

FPB2 

FPB3 

•1 

FPB4 

N 

FPB5 

n 

FPB6 

•• 

MSB 

FPB7 

Not  Used 

FPAO 

Y  Address 

LSB 

FPAl 

II 

FPA2 

ft 

FPA3 

II 

FPA4 

II 

FPA5 

•1 

FPA6 

•1 

MSB 

FPA7 

Not  Used 

SPA7 

WriteNErase  Control 

SPARE 

Not  Used 

GROUND 

Ground 

GROUND 

Ground 

PT 

XY  Pulse 

PC 

Coil  Pulse 

SPARE 

Not  Used 

SPARE 

Not  Used 

SPARE 

Not  Used 

SPARE 

Not  Used 

MSB  :  most  significant  bit. 

LSB  :  least  significant  bit. 

*note:  J5  connector  is  numbered  with  odd 
pins  on  the  outside  and  even  pins 
on  the  inside. 


■  V*^***I#*V*** 
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The  seven  X  addresses  (PPBO-6)  and  seven  Y  addresses 

7 

(FPAO-6)  correspond  to  2  «  128  rows  and  coluBins,  respec¬ 

tively.  For  example,  FPBO-6  >  000  0001  (binary)  corresponds 
to  row  1;  FPAO-6  *  111  1111  (binary)  corresponds  to  column 
127  (rows  and  columns  are  numbered  0-127) . 


Litton  signals 

A  pixel  is  written  when  the  X  and  Y  address  lines  are 
asserted  and  a  positive  logic  level  is  applied  to  the 
write/erase  control  line  (SPA7).  SPA7  latched  high 
(positive  logic  level  -  binary  one)  causes  the  coil  current 


pulse  (PC)  to  go  in  the  clockwise  direction,  which  turns  the 


K" 


nucleated  pixel  on.  If,  on  the  other  hand,  SPA7  is  latched 
low,  the  coil  current  pulses  in  the  counterclockwise  direc¬ 
tion  causing  the  nucleated  pixel  to  turn  off. 

The  timing  of  the  above  mentioned  signals  is  extremely 
critical,  as  shown  in  Figure  2.3.  TPT  and  TPC  are  the 
outputs  of  an  SN74LS123  retriggerable  monostable  multivi¬ 
brator  located  on  the  X  driver  card  and  are  triggered  by  the 
positive  edge  of  inputs  PT  and  PC,  respectively.  A  variable 

resistor  is  connected  between  R.  ^  and  V  on  the  SN74LS123 

xnt  cc 

to  obtain  variable  pulse  widths  for  TPT  and  TPC.  The  pulse 
widths  were  set  as  shown  in  Figure  2.3  to  obtain  optimum 
results.  (Excessive  width  of  TPT  will  cause  LIGHT  MOD  lines 
to  fuse  open  -  the  pulse  is  typically  l//sec  wide)  . 

TPT  must  be  pulsed  while  the  X  and  Y  addresses  and  SPA7 
are  valid  in  order  to  nucleate  a  pixel.  TPC  is  then  pulsed 
to  saturate  the  given  pixel. 
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All  the  required  signals  could  be  obtained  front  the 
Semetex  board  shown  in  Figure  2.4.  The  PT  timing  is 
inherent  in  the  original  hardware  and  the  PC  timing  is 
obtained  through  thesis  software. 

OA00-DA06  were  chosen  to  function  as  the  Y  addresses 
and  DD0-DD6  were  chosen  to  function  as  the  X  addresses;  DD7 
to  generate  SPA7;  WRITE  for  PT;  and  PRINT  for  PC.  The  board 
was  set  to  the  hexadecimal  offset  segment  address  BOOO  by 
properly  setting  the  eight  switches  of  the  interface  board. 
The  high  order  bit  is  on  the  right  and  a  switch  in  the  down 


Pixel  selection  is  accomplished  by  sending  the  desired 
X-Y  address  to  DA00-DA06  and  DD0-DD6.  For  example,  a 
write  to  decimal  50,  100  would  nucleate  the  pixel  at  column 
50,  row  100  (a  write  is  accomplished  in  the  software  by 
poking  data  to  an  address  — *  software  is  explained  in 
Appendix  A) .  DD7  (SPA7  on  the  Litton  X  driver  card)  has  to 
be  high  {binary  one)  to  turn  the  pixel  on.  Before  writing 
to  the  desired  location,  a  logical  OR  of  the  row  (DD0-DD6) 
and  hex  80  (binary  1000  0000)  must  occur.  The  logical  OR 
assures  that  DD7  is  high  without  changing  the  data  in 
DD0-DD6.  If  the  pixel  is  to  be  turned  off,  a  simple  write 
to  column, row  is  sufficient  since  DD7  is  low  (binary  zero) 
by  default. 

Any  write  to  column,  row  triggles  a  WRITE  pulse  since 
WRITE  is  tied  to  DA00-DA06  through  flip-flops  and  the  clock. 

A  write  to  decimal  offset  2049  after  nucleation  of  all 
desired  pixels  will  produce  the  PRINT  pulse,  sending  current 
through  the  coil  and  completing  the  pixel  write/erase 
process.  A  summary  of  pixel  write/erase  operation  is  shown 
in  Figure  2.5.  The  Cl  cable  connection  between  the  Semetex 
interface  board  and  the  Litton  X  driver  card  is  shown  in 


Table  2.2. 


offset  segment 
(Hex  BOOO) 


pixel  on:  logical  OR  of  X  (row)  and  Hex  08 
ixel  off:  X  (row) 


write  X,  Y 
X  -  row 
Y  -  column 


WRITE  IS  pulsed  whenever  Y<255 

-  will  always  be  the  case 
since  0<Y£127. 

-  WRITE  pulse  corresponds 
to  PT  pulse  on  LIGHT 
MOD. 


write  to 
offset  2049 

-  PRINT  pulse  **  PRINT  pulse  corresponds 

- - -  to  PC  (coil  current)  on 

LIGHT  MOD. 


D 


Figure  2.5. 


Pixel  write/erase  control  through 
Semetex  interface  board  to  LIGHT 
MOD  device. 


'V? 


% 


Semetex 

interface  card 

Litton 

X  driver  card 

Signal 

name 

m 

Signal 

name 

7 

DAOO 

1 

FPBO 

19 

DAOl 

2 

FPBl 

20 

DAO  2 

3 

FPB2 

5 

DAO  3 

4 

FPB3 

msm 

DAO  4 

5 

FPB4 

DAOS 

6 

FPB5 

■B 

DAO  6 

7 

FPB6 

22 

DDO 

9 

FPAO 

23 

DDl 

10 

FPAl 

24 

DD2 

11 

FPA2 

12 

DD3 

12 

FPA3 

11 

DD4 

13 

FPA4 

10 

DD5 

14 

FPA5 

9 

DD6 

15 

FPA6 

21 

DD7 

17 

SPA7 

4,25 

GND 

20 

GND 

14 

BUSY 

19 

GND 

8 

WRITE 

21 

PT 

13 

PRINT 

22 

PC 

Cable  signal 
function 


X  Address  Li 


Y  Address 


Write/Erase  Control 

Ground 

Ground 

XY  Pulse 
Coil  Pulse 


Polar-Logarithmic  Coordinate  Transformation 


CGH  for  Scale  and  Rotation  Invariance 

In  this  chapter,  the  scale  and  rotation  invariant 
polar-logarithmic  (lnr-0)  coordinate  transformation  (CT)  is 
used  to  achieve  an  in-plane  distortion  invariant  feature 
space.  The  CT  is  performed  by  a  computer  generated  holo¬ 
gram  developed  on  a  laser  printer  and  a  transforming  lens. 

The  first  section  develops  the  necessary  equations  to 
produce  the  CT  CGH.  The  results  obtained  by  placing  the 
CGH  in  the  input  plane  of  an  optical  processor  are  presented 
and  discussed  in  the  final  part  of  this  chapter.  Appendix  B 
presents  the  fabrication  method  used  to  develop  the  CGH. 

Design  of  the  Coordinate  Transformation  CGH 

The  coordinate  transformed  image  appears  in  plane  of 
Figure  3.1  when  the  CGH  is  placed  directly  against  the 
input  image  in  plane  P^  (6)  .  The  CGH  has  an  ideal  binary 


grating-like  pattern  which  simulates  a  continuous  phase 
transmittance  given  by  h{x,y)  =«exp[  j0  {x,y)  ]  ,  where  0(x,y)  is 
the  phase  distribution.  (The  effect  of  the  phase  transfor¬ 
mation  exists  in  the  first  diffracted  order  of  the  grating) . 


The  Fourier  transform  of  the  product  of  the  input,  f(x,y). 


and  transmittance,  h(x,y) ,  is  performed  by  lens  L^.  In 
plane  this  product  is  defined  mathematically  by 

P(u,v)  «  £  (x,y)exp(j^  (x,y)  ] 

*exp[- j  (2w/\f (xu+yv)dxdy  (1) 

where  i  is  the  wavelength  of  light  used  and  f^  is  the  focal 
length  of  lens  L^. 

By  taking  the  partial  derivative  of  the  Fourier 
kernel,  it  is  found  that  the  coordinate  transformation  which 
maps  an  input  point  (x,y)  to  an  output  point  (u,v)  is 
(7:3100) : 


if  a0(x,y) 

X  - »  “  =  2^  —3^ -  (2) 

if  3^{x,y) 

y  —  *  ?r  Siy) 

To  perform  the  lnr-0  CT,  the  appropriate  expression  for 
0(x,y)  must  be  determined  which  satisfies  (2)  and  (3).  This 
transformation  can  be  expressed  by 

u  *  ln(x^  +  (4) 

V  »  -tan"^(y/x) 


(5) 


and  (5)  in  (3)  then  solving  the  following  equations  for 
0(x,y)  (7:  3099-3104). 


ln(x2^y2)l/2 


-tan~^ (y/x) 


(x,y) 

(6) 

2n 

9  (x) 

3^ (x,y) 

(7) 

2n 

5  (y) 

The  result  is  the  continuous  phase  function 


0(x,y)  “  [2w/\f^]  *  [In  (x^  +  -  ytan  ^(y/x)-x]  (8) 


The  use  of  Equation  (8)  to  form  a  binary  grating  pat¬ 
tern  (CGH)  to  perform  the  lnr-0  coordinate  transformation, 
the  printing  of  that  pattern  on  a  laser  printer,  and  the 
method  used  to  reduce  the  CGH  onto  a  high  resolution  film 
plate  is  discussed  in  Appendix  B. 


Experiment  and  results 

The  final  CGH  used  had  dimensions  of  lOxlOmm  and  a 
carrier  frequency  of  35  lines/mm  .  The  CGH  was  placed  in 
the  input  plane  directly  against  the  input  image  as  shown 
in  Figure  3.2.  Input  images  consisted  of  various  sized 
squares,  circles,  and  letters.  The  laser  source  used  was  a 
25mW  HeNe  (X=6328nm)  and  the  transform  lens  (L^)  had  a  focal 
length  of  300mm.  The  experimental  setup  is  shown  in  Figure 
3.2. 

The  lnr-0  CT  pattern  was  observed  in  plane  ?2  using  a 
Sony  X-38  CCD  (charge  coupled  device)  camera  and  AT&T  frame 
grabber  on  a  Zenith  248  computer  system.  The  Inr-G  CT  is 


collimated  beam 


t 

pinhole 


inputs 


aperture 


25mw  NEC 
HeNe  laser 
(polarized) 


f «300mm 


Sony 

CCD  Camera 


AT&T  frame 
grabber 
in 

Zenith  248 


Figure  3.2.  Scale  and  rotation  invariant  setup, 


present  in  the  first  diffracted  order  in  plane  P2^  The 
angular  separation  of  diffracted  orders  is  dependent  on  the 


carrier  frequency  of  the  grating  pattern  (a)  and  the  wave- 


.  -1 


length  of  incident  light  {\)  by  ®  ■  sin  (aX) .  For  a»35 


lines/mm  and  Xs0.6328^m,  the  angular  separation,  0,  is  1.27 


degrees.  A  diffraction  pattern  which  is  related  to  the 
pixel  size  of  the  CGH  is  also  present.  The  three  dots  above 
the  lnr-9  CT  shown  in  the  experimental  results  corresponds 
to  this  diffraction  pattern. 

The  transformed  patterns  were  observed  to  remain 
unchanged  for  each  class  of  images  (except  for  horizontal 
and  vertical  shifts) .  A  rotation  of  the  input  image  caused 
a  horizontal  (cyclical,  along  the  0  axis)  shift  in  the  CT 
plane  •  This  cyclical  shift  is  proportional  to  the  input 
rotation  (13) .  The  transformed  pattern  was  found  to 
shift  vertically  (along  the  Inr  axis)  for  input  scale 
changes.  These  results  verified  the  use  of  the  CGH  for  the 
desired  Inr-^  CT.  The  results  are  shown  in  Figures  3. 3-3. 8. 

The  experimental  transformations  fade  out  at  small  r 
relative  to  the  theoretical  due  to  the  way  light  is  trans¬ 
ferred  from  the  x-y  plane  to  the  lnr-0  plane.  A  small 
amount  of  light  from  a  central  polar  region  of  the  x-y  plane 
is  transferred  to  the  same  amount  of  area  in  the  lnr-0  plane 
as  a  large  amount  of  light  from  an  outer  polar  region  of  the 


t5  3p-xnac4^^ 
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Conclusion 


This  chapter  developed  the  necessary  equations  to  form 
the  lnr-0  coordinate  transformation  binary  grating  CGH.  The 
results  from  using  the  CGH  in  the  input  plane  of  an  optical 
processor  to  obtain  scale  and  rotation  (but  not  position) 
invariance  was  e-lso  c-resented.  in  order  to  become  position 
invariant,  the  r’^gnitude  of  Fourier  transform  must  be 

used  as  the  input  to  the  CGK.  However,  the  CGH  cannot  be 
placed  in  the  Fourier  transform  plane  since  the  light  is 
not  a  plane  wave  at  that  point.  With  non-parallel  light 
entering  the  CGH,  the  Inr-^  CT  plane  becomes  a  diffracted 
mess  with  input  images  overlayed  in  a  continuous  fashion. 

To  overcome  this,  the  magnitude  squared  of  the  Fourier 
transform  will  be  recorded  using  a  CCD  camera  and  displayed 
on  the  SLM.  The  SLM,  with  parallel  light  passing  through 
it,  will  become  the  input  of  the  CGH.  Another  option  is  to 
take  a  picture  of  the  Fourier  transform  .and  place  the 
developed  slide  into  the  system  to  act  as  the  input  to  the 
CGH.  The  SLM  will  only  give  a  thresholded  binary  image  of 
the  Fourier  transform  whereas  the  picture  will  contain  the 
grey  levels. 

Chapter  two  described  the  modification  and  interfacing 
which  was  necessary  to  incorporate  the  LIGHT  MOD  into  the 
optical  processor  as  a  SLM.  Chapter  4  will  then  present  the 
position,  scale,  and  rotation  invariant  optical  system  using 
slides  since  the  LIGHT  MOD  broke  while  trying  to  perform  the 
experiments . 
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IV .  Position ,  Scale ,  and  Rotation  Invariant 
Feature  Space 

Introduction 

An  ideal  target  recognition  system  should  be  able  to 
detect  the  presence  of  every  occurrence  of  a  particular 
reference  pattern  in  an  input  scene.  However,  six  major 
distortions  in  the  input  pattern  prevent  a  recognition 
system  from  working  properly  (14) : 


1. 

changes  in  target 

position , 

2. 

changes  in  target 

size , 

3. 

in-plane  rotation 

of  the 

target , 

4. 

out-of-plane  rotation  of 

the  target, 

5. 

changes  in  target 

shape , 

6. 

and  changes  in  scene  clutter. 

This  chapter  investigates  the  first  three  distortions. 
Problem  areas  four,  five,  and  six  are  much  more  difficult  to 
compensate  and  were  not  investigated  in  this  study. 

The  use  of  an  optical  correlator  and  matched  spatial 
filter  for  target  recognition  is  advantageous  due  to  its 
ability  to  process  in  parallel.  However,  the  conventional 
correlator  is  unable  to  recognize  scaled  and  rotated 
versions  of  the  reference  object. 

If  two  functions  to  be  correlated  are  identical  except 


for  a  scaling  factor,  the  resulting  correlation  peak  will 
not  be  as  sharp  as  when  the  scaling  factor  is  one.  Like¬ 
wise,  a  rotation  of  one  of  the  functions  will  quickly  de- 


grade  the  correlation  peak  (2:1652;  15). 

Since  an  optical  correlator  was  to  be  used  to  determine 
whether  a  target  in  the  input  scene  matched  a  reference 
target,  the  input  scene  and  template  had  to  be  transformed 
into  a  position,  scale,  and  rotation  invariant  (PSRI) 
feature  space. 

A  PSRI  feature  space  will  be  developed  in  the  first 
section  of  this  chapter.  Three  transformations  were 
performed  on  the  template  and  input  scene  in  order  to  obtain 
a  PSRI  feature  space.  A  Fourier  transform  and  a  Inr-© 
coordinate  transformation  (CT)  was  used  to  change  any  scale 
and  rotation  difference  into  shifts  along  the  feature  space 
axes.  The  next  section  of  the  chapter  shows  the  optical 
system  which  transforms  an  input  image  into  a  PSRI  feature 
space.  The  PSRI  feature  space  can  then  be  used  as  the  input 
to  an  optical  correlator,  as  will  be  shown  in  Chapter  5. 

The  final  section  presents  results  of  the  optical  transfor¬ 
mation  which  develops  the  PSRI  feature  space. 

It  should  be  noted  that  the  optical  PSRI  feature  space 
implemented  in  this  thesis  was  based  on  the  digital  elec¬ 
tronic  algorithm  developed  by  Kobel  and  Martin  on  the  VAX 
computer  system  at  AFIT  (14) .  Many  of  the  scenes  and 
templates  used  in  the  computer  algorithm  were  placed  on  film 
and  used  as  the  inputs  of  the  optical  system.  Comparisons 
to  the  digital  electronic  algorithm  will  be  made  where 
applicable . 


Position  Invariance 

Position  invariance  is  accomplished  by  using  the  mag¬ 
nitude  of  the  Fourier  spectrum,  The  Fourier 

X  y 

transform  (FT)  shift  theorem  will  show  position  invariance 


(10:9)  ; 


f|i(x,y)j-  = 


I(fx'fy)»  then 


F|i(x-a,y-5)|  =  I  (f  f  y )  exp  [-j2n  (af^+  /?fy)j  (9] 


where 


±{x-a,y-l3)  = 


'{  } 


input  image  shifted  by  a-units 
in  the  x-direction  and  5-units 
in  the  y-direction. 

two-dimensional  Fourier 
Transform. 


"This  simply  states  that  a  shift  in  the  input  scene  contri¬ 
butes  a  linear  phase  factor  to  the  Fourier  transform  but 
does  not  impact  the  magnitude  of  the  Fourier  transform  " 
(14:17).  Therefore,  by  using  the  magnitude  (or  magnitude 
squared  of  the  Fourier  spectrum)  there  are  no  variations  due 
to  a  shift  (relocation  of  target)  in  the  input  scene  plane, 
the  magnitude  squared  of  the  FT  was  used  in  the  optical 
system  since  the  intensity  of  the  FT  was  detected. 


Scale  and  Rotation  Invariance 


The  rotation  of  an  object  in  the  input  scene  plane  will 


result  in  an  identical  rotation  of  the  magnitude  of  the 
object's  FT.  The  rotation  of  the  magnitude  spectrum  in  a 


r>Qr. 


rectangular  coordinate  system  can  be  transformed  into  a 
linear  shift  by  mapping  the  magnitude  spectriim  to  a  polar 
coordinate  system  as  shown  in  Chapter  3.  The  magnitude 
spectrum,  |l(f„,f„)  I,  will  now  be  in  the  form  |l(f«,f/j)  I, 

X  y  *  r  “  ' 

with  spatial  frequencies  in  the  angular  (f^)  and  radial  (f  ) 

r 

directions.  The  new  spatial  frequency  coordinates  are  given 
by  (14:17): 


(10) 


(11) 


Any  rotation  will  now  appear  as  a  linear  shift  of  the  polar 
magnitude  spectrum  along  the  angular  axis,  f^. 

A  logarithmic  scaling  of  the  radial  axis  of  the  polar 
magnitude  spectrum  plane  is  then  required  for  scale 
invariance,  "The  Fourier  transform  'similarity*  theorem 
states  that  a  magnification  or  shrin)cage  of  the  coordinates 
in  the  space  domain  corresponds  to  a  respective  5hrin)cage  or 
magnification  of  the  coordinates  in  the  frequency  domain" 
(14:18).  The  similarity  theorem  is  expressed  by  (10:9): 


{ 


F  i  i(x/a,  y/li) 


} 


|a/3|l(af 


(12) 


where  a  and  d  are  the  scaling  constants  for  the  x  and  y 
coordinates . 

The  spatial  frequencies  f  and  f  are  uniformly  scaled 

X  y 

when  (that  is,  the  object  is  uniformly  scaled  in  the  x 
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and  y-direction)  .  As  a  result,  ^  j  "  [  ]' 

Therefore,  the  angular  axis  (£^)  of  the  polar  magnitude 
spectrum  plane  will  not  be  affected  by  a  uniform  scaling  due 


to  the  ratio 


I  1  in  Equation  (10) 

L  J 


by  a 


By  scaling  the  radial  frequency  axis  in  Equation  (11) 


2,2  «2,  2 
a  f  +  13  f 
X  y 


and  talcing  the  natural  logarithm  yields: 


-  (1/2)  In  [  J  a^f^  +  /?^fy^  j 

«  ln(a)  +  (l/2)ln[^  f^  +  f^  j 


The  second  term  of  Equation  (13)  is  just  the  natural 
logarithm  of  the  unsealed  radial  spectrum.  This  shows  that 
a  uniform  scaling  of  an  image  by  a  will  cause  a  linear  shift 
of  the  amount  In (a)  along  the  logarithmic  radial  frequency 


axis . 


An  optical  scale  and  rotation  invariant  feature  space 
was  developed  and  shown  in  Chapter  3  using  the  lnr-0  CT  CGH 
and  a  transforming  lens.  That  same  processor  was  incorpor¬ 
ated  into  the  PSRI  system  by  placing  the  magnitude  squared 
of  the  FT  against  the  Inr-©  CT  CGH. 


Optical  Processor  for  PSRI  Feature  Space 


1L3 


The  experimental  set-up  used  to  transform  an  input 
scene  or  template  from  an  x-y  system  to  a  PSRI  feature  space 
is  shown  in  Figure  4.1. 

The  first  step  in  the  optical  processor  accomplished 

position  invariance.  The  template  or  scene  was  placed  in 

the  input  plane  Pj^  and  Fourier  transformed  by  lens  .  The 

2 

magnitude  squared  of  the  FT  (|FT|  )  was  detected  in  plane 

2 

by  a  Sony  CCD  camera.  The  |FT|  remained  unchanged  when  the 
target  was  relocated  (shifted,  but  not  rotated  or  scaled)  in 
the  input  plane;  thus  giving  the  position  invariance. 

It  was  necessary  to  place  a  high  pass  filter  adjacent 
to  plane  P^  in  order  to  obtain  any  information  from  the 
I  FT  I  .  The  dc  and  lower  order  harmonics  would  cause  the  CCD 
camera  to  wash  out  since  they  were  more  intense  than  the 
higher  frequency  components.  Plus,  there  is  no  real  indivi¬ 
dual  target  characteristic  information  in  the  lower 
harmonics  so  it  was  not  necessary  to  include  them.  Kobel 
and  Martin  (14:31-33)  used  a  bandpass  mapping  function  to 
exclude  undesirable  low  and  high  spatial  frequencies.  They 

established  an  f  .  and  f  by  determining  which  portion 

iDxn  indx 

of  the  spativl  frequency  spectrum  was  generated  by  the  input 

target  silhoette.  They  evaluated  the  size  relationships  of 

the  aircraft  template  silhoette  to  the  overall  template 

window  and  the  size  relationship  between  the  possible  target 

and  scene  window.  The  lower  bound  (f  .  )  was  determined  by 

min 

the  targets  largest  physical  feature  and  the  upper  bound 
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was  determined  by  the  targets  smallest  physical 
feature.  They  found  that  the  fundamental  harmonic  of  the 
largest  physical  feature  (for  a  "typical”  target)  cor¬ 
responded  to  the  4^^  harmonic  whereas  the  fundamental 
harmonic  of  the  smallest  physical  feature  was  the  55^^ 

harmonic.  Therefore,  they  used  a  bandpass  with  f  .  =4  and 

min 

f  =55.  The  results  they  obtained  during  correlation  was 

ItioX 

very  good  so  I  suspected  the  same  would  be  true  using  the 
optical  processor.  A  filter  which  blocked  the  first  4 
harmonics  of  the  |FT|  would  be  the  most  desirable. 

2 

The  program  TARGASLM  was  then  used  to  capture  the  |FT| 

of  the  input  target  and  transfer  it  to  the  SLM  in  plane 

where  the  |FT|  would  be  the  input  to  the  lnr-6  CT  pro- 

cessor.  A  picture  of  the  |FT|  was  also  taken  from  the 

high  resolution  monitor  using  a  35mm  camera  with  200  speed 

color  slide  film.  Since  the  SLM  was  only  able  to  display  a 

binary  sampled  version  of  the  |FT|  ,  the  slide  would  capture 

the  complete  gray  scaled  image  (limited  by  the  resolution  of 

the  CCD  detector  and  film)  but  using  the  slides  would  negate 

the  real-time  capability  of  the  system.  It  was  during  these 

experiments  that  the  SLM  failed  so  all  results  shown  at  the 

end  of  this  chapter  used  the  35mm  slide  technique. 

The  scale  and  rotation  invariant  (lnr-0)  optical 

2 

processor  used  the  |FT|  as  the  input  to  the  lnr-0  CT  CGH. 
Lens  L^  and  a  low  pass  filter  were  used  to  filter  out  the 
harmonics  of  the  SLM  which  were  due  to  the  repeated 
electrode  structure  of  the  SLM  device.  A  very  fine  grid 


pattern  is  formed  by  the  128x128  crossed  electrodes  which 

created  a  large  sine  pattern  in  both  the  x  and  y  directions. 

Only  the  central  harmonic  needed  to  be  passed  since  most  of 

the  displayed  image  is  contained  in  it.  Lens  imaged  the 

I  ft]  onto  the  center  of  the  CGH.  The  CGH,  along  with  lens 

2 

then  performed  the  desired  Inr-d  CT  on  the  |FT|  of  the 
input  target  and  the  resulting  PSRI  feature  space  was 
generated  in  plane  P^.  The  magnitude  squared  of  the  PSRI 
feature  space  was  detected  by  camera  C2  where  it  could  then 
be  transferred  to  the  input  of  an  optical  correlator.  The 
optical  correlator  is  discussed  in  Chapter  5.  A  discussion 
of  the  experiment  and  the  results  are  presented  in  the  next 
section . 

Experimental  results 

The  goal  of  these  experiments  was  to  transform  an  input 
template  or  scene  into  a  PSRI  feature  space.  With  such  a 
featue  space,  a  target  shifted  to  another  part  of  the  scene 
would  have  the  same  shape  after  transformation  as  the 
template.  Any  rotation  or  scale  change  of  the  target  with 
respect  to  the  template  would  relate  to  a  simple  horizontal 
or  vertical  shift  in  the  PSRI  feature  space.  Templates 
consisted  of  geometrical  shapes,  letters,  and  aircraft 
silhoettes.  The  aircraft  silhoettes  were  those  used  by 
Kobel  and  Martin  in  their  digital  electronic  computer 
algorithm.  They  were  obtained  from  video  files  stored  on 
the  VAX  computer  system  and  placed  on  200  speed  35mm  slides 


1 
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to  be  used  as  inputs  to  the  optical  system.  The  real  scenes 
containing  the  airplanes  were  also  obtained  from  the 
computer  system. 

The  first  set  of  experiments  consisted  of  placing 
simple  geometric  objects  in  the  input  plane  and  per¬ 
forming  the  FT  with  lens  By  using  the  high  pass  filter 

in  front  of  camera  quality  FT*s  were  observed  in  plane 
P2.  TARGASLM  was  then  used  to  capture  the  |FT|  in  video 
memory  and  transfer  the  |FT|  to  the  SLM  in  plane  P^  to 
function  as  the  input  to  the  CGH.  Best  results  were 
obtained  when  only  the  center  256x256  of  the  512x400  video 
image  was  sent  to  the  SLM.  The  threshold  value  remained  as 
a  variable  input  to  the  user  since  the  amount  of  energy 
passing  through  the  input  varied  greatly.  At  the  time  an 
8mW  laser  was  used  as  the  source  and  the  CCD  camera  could 
barely  detect  the  PSRI  feature  space  in  the  diffracted  order 
of  the  CGH/lens  output  (P^).  The  system  appeared  to  work 
properly  but  the  SLM  failed  irremedially  before  any  useful 
results  were  obtained  after  incorporating  the  25mH  laser 
into  the  system  as  the  source. 

The  first  set  of  results  consisted  of  a  small  square 

and  a  larger,  rotated  square  as  the  inputs  (Figure  4.2). 

2 

The  I  FT  I  of  these  two  inputs  are  shown  in  Figure  4.3.  The 
2 

I  FT  I  were  also  placed  on  slides  and  inserted  into  plane  P^ 
where  they  were  imaged  onto  the  center  of  the  lnr-0  CGH. 

The  resulting  PSRI  feature  space  detected  in  plane  P^  is 
shown  in  Figure  4.4.  Note  the  vertical  (along  the  Inr  axis) 


and  horizontal  (0  axis)  shifts  associated  with  the  scale  and 
rotation  difference  between  the  two  inputs.  With  Figure 
4.4a  used  as  the  reference  template  (matched  spatial  filter) 
and  Figure  4.4b  used  as  the  input  to  an  optical  correlator, 
it  is  presumed  that  a  very  sharp  correlation  pea)c  would 
result  after  a  crosscorrelation  is  performed. 

In  the  next  set  of  results  presented  in  this  Chapter, 
the  side  view  silhoette  of  an  F15  was  used  as  the  template 
and  a  scaled  and  rotated  version  of  the  F15  was  used  as  a 
comparison  (Figure  4.5).  Also  included  with  the  analysis 
was  an  F15  amongst  some  clouds  in  a  real  scene  (Figure  4.6). 

The  PSRI  feature  spaces  for  both  the  templates  and  the 
real  scene  are  shown  in  Figure  4.7.  Again,  this  feature 
space  shown  could  be  used  as  an  input  to  an  optical  cor¬ 
relator  since  each  only  differs  by  a  horizontal  or  vertical 
shift . 

When  the  airplane  templates  and  real  scenes  were  used, 
the  results  were  not  as  clear  as  those  obtained  using  simple 
geometrical  objects.  This  can  be  attributed  to  three 
observed  problem  areas.  The  first,  and  most  obvious,  is 
that  the  objects  are  more  complicated  malcing  the  results 
less  intuitive  (there  is  really  no  a  priori  )tnowledge  of  how 
the  PSRI  feature  space  should  loo]()  . 
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The  next  problem  concerned  the  detection  of  the  |FT| 
after  passing  through  the  high  pass  filter.  The  filter 
consisted  of  an  in)c  blot  in  the  middle  of  a  glass  slide. 


Figure  4.6. 


F15  scene. 


The  ink  blot  did  not  act  as  a  pure  block  of  the  first  few 
harmonics  but  seemed  to  block  in  some  areas  and  only 
attenuate  others.  More  importantly,  the  filter  was  non- 
circular  which  caused  the  [FT's]  to  be  non-symmetric .  A 
very  serious  problem  arose  since  a  characteristic  of  the  FT 
is  there  inherent  symmetry. 
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The  last  and  most  troublesome  problem  was  in  the 

centering  of  the  |FT|  slides  onto  the  CGH.  For  the  Inr-d 

CT  to  exist  in  the  diffracted  order  of  plane  the  input 

object  absolutely  must  be  centered.  Since  each  |FT|  was 

mounted  slightly  different  in  its  respective  slide,  they 

all  had  to  be  individually  centered  on  the  CGH.  That  was  a 

very  subjective  process  when  the  PSRI  feature  space  result 

was  unknown.  That  extra  variable  led  to  huge  errors  when 

off  center  by  only  a  few  tens  of  microns.  This  would  not 

2 

have  been  a  problem  if  the  SLM  was  used  since  each  |FT| 
would  have  been  detected  in  the  same  location  and  then 
transferred  to  the  SLM  and  displayed  in  precisely  the  same 
location  each  time.  Therefore,  once  the  system  was  centered 
it  would  not  have  to  be  adjusted  again. 

Conclusion 

This  chapter  developed  the  PSRI  feature  space  using  a 
FT  and  a  Inr-G  CT  CGH  in  an  optical  processor.  The  system 
was  shown  to  work  in  real  time  using  a  magneto-optic  spatial 
light  modulator.  The  next  chapter  will  give  background 
theory  on  an  optical  correlator  using  a  matched  spatial 
filter.  Also  included  in  the  next  chapter  will  be  a 
description  of  how  the  PSRI  feature  space  could  be  used  as 
the  input  to  the  optical  correlator  built  for  this  thesis 
effort. 


V. 


Optical  Correlation 


Introduction 

Optical  correlation  has  been  mentioned  throughout  this 
thesis  as  a  method  to  “search"  an  input  scene  for  a  parti¬ 
cular  reference  target.  The  problems  with  the  correlator 
which  arise  from  scale  and  rotation  changes  in  the  input 
with  respect  to  the  reference  (template)  were  used  as  the 
basis  to  transform  the  input  scene  into  a  PSRI  feature 
space . 

The  PSRI  feature  space  of  the  reference  target 
(template)  would  be  used  to  construct  a  matched  spatial 
filter  (MSF)  and  placed  into  the  correlation  plane  of  the 
optical  correlator.  The  PSRI  feature  space  of  a  scene  would 
then  be  used  as  the  input  to  the  optical  correlator.  The 
shape  and  intensity  of  the  resultant  correlation  pea)c 
(spike)  would  disclose  whether  the  reference  target  was 
located  in  the  scene. 

Although  no  MSF's  were  constructed  and  no  correlation 
experiments  performed,  the  inclusion  of  this  chapter  is  for 
completeness  of  the  optical  pattern  recognition  system.  The 
first  section  presents  the  theory  of  MSF's  and  optical 
correlators  using  the  method  of  A.  Vander  Lugt  (10:171-192). 

The  next  section  explains  the  experimental  configuration 
that  was  constructed  in  this  thesis  effort.  Again,  no 
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experimental  results  using  the  optical  correlator  were 
obtained  due  to  a  time  limitation. 

Theory 

A  simplified  optical  correlator  is  shown  in  Figure  5.1. 
An  input,  g(x^,yj^),  is  placed  in  the  front  focal  plane  i?^) 
of  lens  L^.  Lens  takes  the  FT  of  the  input  producing 
G(f  ,f  )  at  plane  P-.  A  mask  of  transmittance  H(f^,f,  )  is 
placed  in  plane  P2  to  function  as  the  MSF  (10:178).  The 
product  of  G  and  H  exists  immediately  after  plane  P2  and 
lens  L2  takes  the  FT  of  this  product  and  produces  F{GH|  in 
plane  P^. 

The  correlation  between  two  functions,  g  and  h,  can  be 

* 

found  by  taking  the  FT  of  the  product  GH  where  G=F(gl  and 
H=F{hl  (10:175).  By  creating  a  mask  of  intensity  H* 

(complex  conjugate  of  H)  in  plane  P2  of  Figure  5.1,  the 
optical  processor  is  performing  the  correlation  between  the 
input  g  and  a  reference  function  h.  If  the  input  function 
is  identical  to  the  reference  function,  then  the  output  will 


input 

(Sf (Xj^,yj^) ) 


correlation 

plane 


Figure  5.1.  Optical  correlator  using  MSF 


SB? 

be  the  auto  correlation,  and  an  intense  spike  will  appear  in 
plane  .  Otherwise,  a  cross  correlation  will  be  performed 
resulting  in  a  spike  in  plane  whose  intensity  is  depen¬ 
dent  on  the  similarity  between  the  input  and  reference. 

The  MSP  can  be  created  using  a  Vender  Lugt  filter  which 

< 

incorporates  an  interferometric  system  as  shown  in  Figure 
5.2.  The  desired  impulse  response,  reference  function  (h) . 
is  placed  in  plane  P^  and  Fourier  transformed  by  lens  L^. 

The  FT  of  h  is  interfered  with  a  plans  wave  incident  on  the 
film  in  plane  P.  at  an  angle  9. 


Figure  5.2.  MSF  construction  using 

Vander  Lugt  filter  (10:172)  . 


The  total  intensity  distribution  on  the  film  is 
determined  by  the  two  amplitude  distributions  present, 
yielding  (10:172): 


KXj.yj)  -  |r^.xp(-32«ayj)  +  | 

■  ]|  *  if  ]exp(j2nay2) 


*r*2  ^21 

Tf  «  [If 'XfJ®^P<-32"ay. 


where  a  =  is  the  amplitude  of  the  plane  wave,  and 

f  is  the  focal  length  of  lens  The  exposed  film  is  then 

developed  to  produce  a  transparency  with  an  amplitude 

transmittance  proportional  to  the  intensity  distribution 
incident  during  exposure  (10:172): 

t(X2,y2)  «  +  ^|H|^  +  if  H  exp(j2nay2) 

^o  * 

+  H  exp(-32nay^)  (15) 

The  film  plate  is  then  placed  into  plane  ^2  optical 

correlator  (Figure  5.1)  to  function  as  the  MSP.  With  g(x,y) 
as  the  input,  the  field  directly  behind  the  mask  is  (10:175) 


A  I  A  I 

^O 

+  —5—0  H  G  exp(-j2way  ) 

Lens  of  Figure  5.1  performs  the  FT  on  the  field 
described  in  Equation  (16)  and  produces  (10:175) 


2  1  * 

^  <  "X^ , -y3) '3  (x^  ,  y3 )  ] 

^  £ 
r 

+  yI (X3'y3) *«<^3'y3>  ix^,y^+a\t) ] 

■*■  (-X3,-y3)  *3(X3,y3)  *6  (X3,y3-a\f )  ]  (17) 


in  plane  P3.  The  fourth  term  of  this  equation  is  the  cross 
correlation  of  q  and  h  centered  at  (0,a\f)  and  the  one  of 
concern  in  an  optical  correlator.  Certain  procedures  must 


be  followed  when  constructing  the  filter  to  assure  that  the 
other  terms  in  the  output  do  not  overlap  with  the 
correlation  term  (10:175). 


Experimental  Configuration  for  Optical  Correlator 

The  experimental  configuration  for  optical  correlation 
using  a  MSF  is  shown  in  Figure  5.3.  The  correlator  was 
built  into  the  existing  PSRI  transformation  processor  by 
adding  a  beamsplitter  (BS^)  to  create  a  reference  beam.  The 
PSRI  feature  space  would  be  place  in  the  input  plane  P^  and 
a  MSF  would  be  either  created  or  placed  in  plane  P^  by 
interfering  the  FT  of  the  PSRI  feature  space  with  the 


reference  beam.  To  perform  an  optical  correlation,  the  ex¬ 
posed  MSF  would  be  placed  in  plane  P^  and  the  reference  beam 
would  be  bloclted.  The  correlation  of  the  input  PSRI  feature 


space  with  the  MSF  (of  a  template)  would  appear  in  plane  Pg 


A  similar  configuration  was  used  by  David  Neidig  (15) 
to  make  MSF ’ s  and  perform  optical  correlations.  MSF ' s  could 
be  created  using  the  configuration  in  Figure  5.3  by  placing 
a  high  resolution  film  plate  in  plane  and  interfering  the 
reference  beam  with  the  FT  of  the  input.  However,  before 
any  optical  correlation  could  be  performed,  an  energy 
normalization  technique  would  have  to  be  incorporated  into 
the  system  as  described  by  Neidig  (15:12-14). 

Additional  comments  about  the  optical  correlator  are 
presented  in  the  Recommendations  section  of  the  next 


chapter . 


Conclusions 


Conclusions  and  Recommendations 


The  Litton  Iron  Garnet  H-Triggered  Magneto-Optic  Device 
(LIGHT-MOD)  was  successfully  modified  and  integrated  to 
function  as  a  spatial  light  modulator  (SLM)  in  a  real-time 
optical  processor.  An  environment  was  set-up  which 
integrated  the  LIGHT-MOD  with  a  Zenith  248  microcomputer  and 
AT&T  frame  grabber. 

Secondly,  a  real-time  rotation  and  scale  (but  not 
position)  invariant  optical  system  utilizing  a  computer 
generated  hologram  (CGH)  was  constructed.  The  system  demon¬ 
strated  excellent  performance  against  a  broad  range  of  non- 
cluttered  objects,  targets,  and  scenes. 

A  position,  scale,  and  rotation  invariant  (PSRI) 
optical  system  was  tested  in  the  final  section.  Results 
from  this  system  were  inconclusive  since  correlation  was  not 

performed.  The  magnitude  squared  of  the  Fourier  transforms 

2 

(IftI  )  were  slightly  distorted  due  to  an  imperfect  high 
pass  filter.  After  the  LIGHT-MOD  physically  failed,  the 
centering  of  inputs  onto  the  CGH  became  a  serious  problem 
and  added  distortions/errors  to  the  final  PSRI  feature 
space . 

An  optical  correlator  was  constructed  but  never  tested. 


Recommendations 


For  the  optical  system  to  function  as  a  real-time  PSRI 
pattern  recognition  system  certain  areas  must  be  further 
investigated . 

1)  A  new  magneto-optic  substrate  should  be  acquired 

for  the  LIGHT-MOD.  The  LIGHT-MOD  should  then  be  used  as  a 

2 

SLM  to  display  the  |FT|  of  an  input  onto  the  CGH.  Doing  so 
would  circumvent  the  previously  mentioned  centering  problem 
and  give  a  real-time  capability. 

2)  A  set  of  various  sized  high  pass  filters  should  be 

2 

created  to  improve  the  |FT|  and  get  rid  of  unwanted  har¬ 
monics.  They  can  be  created  using  the  Dekacon  optical 
system  described  in  Appendix  B  by  imaging  circles  on  high 
resolution  film  plates. 

3)  Another  SLM  or  a  liquid  crystal  display  TV  could  be 
used  to  function  as  a  variable  high  pass  filter. 

4)  The  programs  written  in  C  to  integrate  the  LIGHT- 
MOD/Zenith  248/frame  grabber  system  should  be  added  to  and 
improved  upon  to  incorporate  more  functions  of  the  AT&T 
frame  grabber,  increase  the  display  speed  of  the  LIGHT-MOD 
when  used  in  conjunction  with  the  CCD  cameras  and  frame 
grabber,  and  add  more  functions/options  to  create  a  more 
versatile  system. 

5)  The  optical  correlator  should  be  fully  tested  and 
then  used  to  correlate  the  PSRI  feature  spaces.  An  energy 
normalization  technique  should  also  be  incorporated. 


6)  The  PSRI  feature  spaces  obtained  using  the  optical 
system  should  be  transferred  to  the  VAX  and  correlated  using 
the  Kobel-Martin  algorithm  (14),  Likewise,  the  PSRI  feature 
spaces  obtained  with  the  Kobel-Martin  algorithm  could  be 
used  as  inputs  to  the  optical  correlator. 


Software  Tools 


This  Appendix  describes  and  presents  the  software 
programs  used  in  conjunction  with  the  thesis  effort.  The 
first  section  describes  the  vendor  provided  software:  AT&T 
TARGA  software  which  controls  the  frame  grabber  installed  on 
the  Zenith  248.  Only  brief  descriptions  are  provided  since 
a  full  description  can  be  found  in  the  manual  (1) .  An 
explanation  of  the  required  software  to  display  images  on 
the  Litton  spatial  light  modulator  (SLM)  is  presented  in  the 
second  section.  The  last  section  gives  a  detailed  descrip¬ 
tion  of  the  software  written  to  interface  the  AT&T  frame 
grabber  with  the  Litton  SLM. 

All  software  was  written  in  C  using  the  Computer 
Innovations  C86  compiler.  In  order  to  modify  or  add  to  the 
existing  software,  one  must  become  familiar  with  the  C 
programming  language,  a  C  compiler  (Computer  Innovations  C86 
is  installed  on  the  Zenith  248  which  is  interfaced  with  the 
optical  system) ,  and  a  linker.  The  programs  developed  by 
the  researcher  and  discussed  here  are  included  at  the  end  of 
this  Appendix. 

AT&T  TARGA  System 

The  AT&T  TARGA  board  is  an  analog-to-digital  converter 
for  a  frame  of  video  data  which  is  input  as  composite  video 
from  a  charge  coupled  device  (CCD)  camera.  The  TARGA  system 


then  digitizes  the  data  on  a  0  to  255  grey  scale  level  in  a 
512x400  pixel  array.  The  image  has  a  higher  resolution  at 
the  TARGA  board  since  the  Sony  CCD  array  is  491x384.  These 
digitized  frames  can  be  saved  to  disk  by  using  the  supplied 
AT&T  TARGA  software. 

The  TARGA  board  is  accompanied  by  two  diskettes  of 
software:  (1)  TARGA  software  tools  and  (2)  the  TARGA 
demonstration  disk.  The  software  tools  are  a  library  of  C 
programs  that  utilize  the  TARGA  board  functions.  Each  of 
these  programs  (subroutines)  can  be  compiled  individually  to 
be  called  by  other  programs  or  used  alone.  TESTTARGA  and 
TRUE_ART  are  the  two  demonstration  programs.  TESTTARGA 
allows  the  user  to  run  individual  subroutines  from  the  tools 
library  to  better  understand,  modify,  or  debug  them. 

TRUE_ART  is  a  user  friendly  program  that  utilizes  all  of  the 
TARGA  tools  and  allows  the  user  to  fully  operate  the  camera 
and  frame  grabber.  TRUE_ART  is  a  very  versatile  program 
which  was  used  greatly  during  the  thesis  effort  to  view 
optical  results,  capture  and  store  frames,  write  text  on 
images,  and  draw  on  images. 

All  of  the  TARGA  programs  reside  on  the  hard  disk  in 
the  directory  Targa  and  TARGA  software  tools  are  in  the  sub¬ 
directory  TargaNTools.  The  program  TRUB_ART.BXE  is  invoked 
by  running  TARGA8.BAT  in  the  Targa  directory.  An  AT&T  logo 
appears  on  the  screen  and  the  user  is  instructed  to  press 
the  Enter  key.  A  menu  is  displayed  and  selections  are  made 


by  moving  the  cursor  with  the  up/down  cursor  keys  to  the 
desired  function  and  pressing  the  FIO  key.  FIO  and  F9  are 
used  to  invoke  all  the  TARGA  commands,  not  the  Enter  key. 

FIO  is  used  to  execute  a  function  and  F9  is  pressed  to  go 
back  one  menu.  It  should  be  noted  that  the  cursor  can  be 
made  to  move  through  the  menus  faster  by  pressing  the  Scroll 
Lock  key.  See  the  TARGA  user's  manual  for  a  description  and 
operation  of  all  commands  and  their  functions  (1) . 

LIGHT  MOD  Software 

The  first  program  written  for  the  LIGHT  MOD  was 
RASTER. C.  This  program  erases  the  LIGHT  MOD  display,  writes 
horizontally  to  the  display  (15  rows  at  a  time),  erases  the 
display,  writes  vertically  to  the  display  (15  columns  at  a 
time) ,  then  erases  the  display  again.  This  sequence  is 
followed  20  times  before  the  program  is  exited.  RASTER. C 
was  used  to  test  and  fine  tune  the  LIGHT  MOD  after  the 
interface  and  modification  was  performed.  While  this 
program  was  running,  the  X  and  Y  current  levels,  print 
pulse,  and  coil  pulse  were  adjusted  until  each  pixel  on  the 
display  was  turned  all  the  way  "on"  and  "off",  thus 
obtaining  the  greatest  dynamic  range.  A  full  frame  could  be 
displayed  in  approximately  300msec. 

The  full  display  is  erased  by  sending  an  X-Y  current 
pulse  to  every  pixel  followed  by  a  counterclockwise  coil 
pulse.  A  pixel  is  accessed  (turned  on  or  off)  in  the 


software  by  poking  data  to  a  specific  offset  address.  The 
data  represents  a  row  on  the  SLM  and  the  offset  address 
represents  a  colxisin  (row  and  column  both  can  take  on  values 
from  0  to  127).  The  Write  pulse  is  activated  any  time  the 
offset  address  is  less  than  2049.  Therefore,  anytime  a 
column  is  addressed  the  Write  pulse  is  generated  to  select 
that  pixel  and  place  it  in  its  neutral  state.  Any  poke  to 
offset  address  2049  will  cause  the  coil  to  pulse  which 
saturates  any  pixel  in  the  neutral  state  and  completes  the 
write  or  erase  function. 

In  C  language,  the  command  to  poke  data  to  certain 
addresses  is 

pokeb (of fset , segment , byte) 

where:  offset  =  the  offset  segment  relative 

to  the  segment  value, 
segment  =  segment  address  of  the  device 
(for  the  LIGHT  MOD  it  is  at 
hex  BOOO) . 

byte  value  or  data  to  be  sent  to 
the  offset. 

The  modification  and  interface  was  built  in  such  a  way  that 
the  offset  value  directly  corresponded  to  a  column  and  the 
byte  value  directly  corresoponded  to  a  row.  The  row  had  to 
be  logically  ORed  with  hex  80  when  a  pixel  was  to  be 
written.  This  logical  OR  operation  sent  SPA7  high  and 
caused  the  next  coil  pulse  to  be  in  the  clockwise  direction. 
When  a  pixel  was  to  be  erased,  the  actual  value  of  row 
(0-127)  was  all  that  was  needed  since  that  defaulted  SPA7  to 
a  low.  With  SPA7  low,  the  coil  pulse  would  be  counter- 


clockwise  causing  the  pixel  to  turn  off.  The  program 
RASTER. C  shows  exactly  how  to  write  or  erase  each  and  every 
pixel  in  the  LIGHT  MOD  display.  Any  other  access  to  the 
LIGHT  MOD  display  is  a  variation  of  this  routine. 


Interface  Software 


Once  the  LIGHT  MOD  was  shown  to  be  operational,  it  was 
necessary  to  interface  the  device  with  the  Zenith  248  in 
order  to  display  real  scenes.  The  remainder  of  the  programs 
mentioned  in  this  Appendix  allow  one  to  manipulate  the 
camera/frame  grabber  system  and  display  stored  or  live 
images  on  the  LIGHT  MOD. 

TARGASLM.C  is  the  user  friendly  program  which  inter¬ 
faces  the  AT&T  TARGA  system  frame  grabber  to  the  LIGHT  MOD. 
This  program  initializes  the  system  and  presents  self 
explanatory  menus  which  can  be  chosen  to  operate  the 
frame  grabber  and/or  display  images  on  the  LIGHT  MOD. 
TARGASLM.EXE  is  located  on  drive  C  (the  hard  disk)  under  the 


directory  Targa  on  the  Zenith  248  which  is  interfaced  with 
the  optical  system.  This  and  all  other  programs  developed 
for  this  thesis  are  also  on  a  5  1/4  "  floppy  diskette 


labeled  SLM_PROGRAMS . 

The  desired  function  in  TARGASLM  is  executed  by  typing 
the  associated  letter  (shown  in  the  left  column  of  the  menu) 


followed  by  the  enter  key. 

This  program  calls  various  AT&T  TARGA  programs  from  the 


TARGA  library  as  well  as  two  programs  to  display  images  on 
the  LIGHT  MOO:  CCDSLM.C  and  SLM.C.  All  TARGA  programs  are 
described  in  the  AT&T  TARGA  manual  (1) . 

CCDSLM.C  and  SLM.C  are  very  similar  programs  which  are 
derived  from  the  TARGA  program  GETPIC.  In  these,  an  image 
file  from  disk  is  opened,  converted  to  a  128x128  binary 
image  file,  and  sent  to  the  LIGHT  MOD.  Both  have  a  variable 
thresholding  value  which  is  input  by  the  user  to  convert  the 
255  grey  scale  TARGA  image  to  a  binary  image  (each  pixel 
either  on  or  off).  The  threshold  value  (0-255)  is  chosen  to 
obtain  the  best  possible  binary  image  on  the  LIGHT  MOD. 
CCDSLM.C  takes  the  full  512x400  TARGA  image  and  sends  every 
4^^  column  and  ever  3*^^  row  to  the  128x128  LIGHT  MOD.  This 
allows  the  full  image  to  be  displayed  on  the  LIGHT  MOD  but 
the  image  becomes  distorted  and  looses  image  resolution  due 
to  the  sampling  technique. 

SLM.C  only  takes  the  center  256x256  of  the  TARGA  image 
and  sends  every  other  pixel  to  the  LIGHT  MOD  in  order  to 
obtain  higher  resolution  of  the  image  on  the  LIGHT  MOD  and 
prevent  distortion  of  the  image. 

All  programs  are  commented  well  and  therefore  only 
described  breifly  here.  Any  modifications  or  additions  to 
TARGASLM  or  the  other  programs  must  be  followed  by  proper 
compilation  and  linking.  It  is  very  important  that  all 
called  programs  are  included  during  the  linking  process. 
After  compiling  the  program  (run  CC.BAT  when  using  C86) ,  the 


r: 
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linker  is  invoked  by  running  LINK.  The  compiled  programs 
object  file  (.OBJ)  must  be  included  followed  by  the  .OBJ 
files  of  all  called  programs  that  are  not  in  a  library.  The 
libraries  C86S2S  and  TARGA  must  also  be  included  since 
C86S2S  contains  all  the  C  programming  functions  and  the 
TARGA  library  contains  the  .OBJ  files  of  all  the  TARGA  tools 
programs.  The  system  compiler  and  linker  are  described  in 
the  appropriate  manuals. 


Progr&a:  RASTER. C 


Written  by;  Michael  w.  Mayo 
20  July  1987 

*  This  prograa  is  used  to  test  the  operation  of  the 

*  Litton  spatial  light  modulator  ISLM).  Rows  and  columns 

*  are  alternately  turned  on  (15  at  a  time)  while  the  dispiav 

*  IS  erased  after  each  horizontal  (rows)  or  vertical  I  cols  i 

*  write. 

*  While  the  program  is  running,  the  polarizer  whicn  is 

<  locatea  after  the  SLM  can  be  rotated  until  the  desired 

*  dynami;  range  is  reached  (that  is.  an  ’on’  pixel  is  as 

*  bright  as  it  can  be  and  an  ‘off  pixel  is  as  dark  as  it 

*  can  be  Other  adjustments  can  be  made  to  the  and 

*  V-current  levels.  WRITE  pulse  width,  and  COIL  pulse  width. 

*  However,  one  must  be  very  careful  when  adjusting  these 

<  parameters.  They  are  set  at  what  I  feel  are  optimum  setti.ngs 

*  and  their  adjustment  should  be  avoided  if  at  all  possible 

*  (refer  to  the  Litton  LIGHT  MOD  manual), 
t  / 

•define  0S_'1ASK  0x80  /*  Logical  OR  with  row  sets  SPAT  high 

«  to  allow  a  pixel  write 

*/ 

extern  unsigned  char  pokeb( ) ; 

unsigned  3egment  =  0xb000 ;  /*  Segment  addrw.>s  of  LIGHT  .MOD  device.  */ 

unsigned  row,  col;  /*  Row  and  Col  of  the  LIGHT  MOD  display.  */ 

char  byte;  /*  Data  sent  to  address.  »/ 

main ;  ) 

( 

1  n  t  i  o  o  c  ; 
int  a  . 

»  Gees  througn  program  20  times.  */ 
for  iloop=0;  loopcDO;  loop'*"^ )  ( 

erase)  /*  Calls  function  to  erase  the  entire  display.  */ 

/* 

«  Section  writes  horizontally  to  the  LIGHT  MOD 

*  15  rows  at  a  time. 

*  / 

printf '  ■ \n\ tBeginning  horizontal  write  to  SLM.\n"): 
a  =  0 ; 

while  a<  128  )  ( 

for'rowsa;  row<=a+l5;  row-t-f  i  ( 
for(col  =  0:  col<128:  col-^-^ii 

byte  =  row  ON_MASK;  /•  Puts  SPA"  high  for 

write  operation  */ 

pokeb(ooi .segment, byte  I ;  /*  PT  pulse  to  nucleate  pixels  «/ 

1 

1 

pokeb(  2049, segment, 253 ) ;  /*  PC  pulse  to  saturate  pixels  «/ 

ay=l6 ; 

I 

printf (  "\n\tHorizontal  write  complete . \n’ ) ; 


•  rase(  ) ; 


/t 

«  Section  write*  vertically  to  the  LIGHT  MOD 

•  15  columns  at  a  time. 

*/ 

•  printf I "XnN tBe*inning  vertical  write  to  3LM.\n  ; 
a  =  0 ; 

while  (  a< 1 28 ) f 

for!rou=0;  row<127:  row«-+)( 

foricolsa:  col<=a*15:  coi*^H 
byte  =  row  !  OS_MASK; 

pokeb( col , segaent . byte ) ;  /*  PT  pulse  to  nucleate  pixels  */ 

1 

\ 

1 

pokebi 2049 , segment , 2S5 ) :  /*  PC  pulse  to  saturate  pixels  */ 

a'«'=  16 ; 

) 

printf ( "\n\tVertical  write  complete . \n" ) ; 

) 

printf ( "\n\nProgram  over!"); 


erase( )  /*  Function  erases  the  entire  SLM  */ 

printf ( "\n\n\tErasing  SLM  display . \n" ) ; 
forlrowsO;  row<sl2T;  row^-f  i  ( 
for(col=0:  col<sl2T;  col*^i  1 
byte  s  row; 

pokeb( col , segment . byte  I ;  /»  PT  pulse  to  nucleate  pixels  »/ 


pOKeoi 2049 . segment , 25  3  .  ;  /»  PC  pulse  to  saturate  pixels  */ 

pr;ntf (  AnVtErase  operation  complete . \n\n  ’ )  : 


/ 


Prfa:  TAACASLM.C  Written  by  Micheel  W.  Mayo.  20  Aug  198T 
Subroutines;  ccdala.  sla 

TAJiCA  subroutines:  Graphlnit,  GetPic.  SetLiveMode . 

GrabFraae.  PutPic,  GraphEnd. 


This  IS  a  user  friendly  prograa  ll  hope :  I  which 
interfaces  the  AT4T  TARGA  f raae  grabber  system  with  the 
Litton  spatial  light  aodulator.  Proa  this  program  you 
can  operate  the  caaera/frame  grabber  and/or  display 
laages  cn  the  spatial  light  sodulator  (SLM).  More 
modifications  could  be  addeo  to  this  program  so  all 
the  TARCa  software  functions  can  be  incorporated.  Oue 
tc  a  lack  of  time  and  necessity,  only  the  above  mentioned 
TARGA  programs  were  incorporated.  To  get  full  use  of  tne 
TaRCA  system,  the  ATIT  provided  program  TRUE_ART.£XE  snuuld 
be  used. 


/ 


1 


I 


•include  "tardev.h" 

•include  <stdio.h> 

extern  struct  TARstruct  *targa; 

•  define  CLS  '’\033[2J"  /*  Escape  code  for  clear  screen  */ 

•define  CUP  ' \033 [ Xd ; XdH'  /»  Moves  cursor  to  position  row.  coi  */ 

main  (  ) 

( 

setup  I  !  : 
menu  \  <  ; 

print:  1  AOSSfOm  /«  ill  attributes  off:  black  on  u-nite 

prir.tf  CLS  . 

I 

setup  1 

i 

printf  (  "XOSO  i' -iA  :  30a  '  i  ;  /  *  Sets  a  blue  oaCKground  with  blacit  letters  * 
pr intf I CLS  j :  / » 

*  Clears  screen  and  positions  cursor 

*  in  upper  left  corner. 

»/ 

pr 1 n t  f 1 CC  P . 7 . 20  1  : 

orintfi  ■  AT4T  FRAME  GRABBER  and  ’  )  :  pr int f ( CUP . 8 . 20  )  : 

printf(  SPATIAL  LIGHT  MODULATOR  OPERATING  SYSTEM  '  . 
pr 1 nt f I  CUP . 2  1  . 30  I  ;  printf (  "Press  ENTER  to  oegin  )  ; 
printfi  "SOSOtSa'  i  :printfiCtP.  16.28)  ;printfi  "NOSOf  36m''  i  ; 
printf ( "Turn  on  High  res  monitor  "  ) ; pr intf ( CUP .  1 7 , 28  I  ; 
printf! "  and  camera  "); 

printf  <  CUP ,23.1); 

i.n  i  le )  ge  tchar  ;  )  '=  '\n');  /*  waiting  for  return  ke>’  */ 


printf! CLS  )  : 

printf  (  CUP.  11.28)  ;  printf  ("  INITIALIZI.NG  TARCA  SYSTEMVn"  )  ; 
if  ( Graph  Ini t ! - I ) :  =  -  i  )  I 

printf!  ''.tNtitin  RASTER  COUNTER  IS  BROKE.NAn" )  ; 
printf!  NtStYt  I/O  SPACE  PROBABLY  WRONC\n\n  "  1 ; 
j>rintf(  'YtYtSt  GRAPHINIT  FAILED  ); 

printf ! "\033(0a  "  )  ;  /»  All  screen  attributes  off  */ 
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■enui  I 


/*  Sets  up  screen  for  asin.Mnu  */ 


t 

printf ( '\033 [ 4 la" ) : printf ( CLS ) :  /*  Red  background;  clear  screen  */ 

printf  (  "XOSS [ 44  ;  32b'' )  ;  /*  Black  window,  green  letters  */ 


printf(CL’P,  3,18);  printf;  "  MAIN  MENU  i 

printf(CUP.  4,18);  printf ( 

pnntflCUP.  3.18);  printf (  ker  function  "i 

printf (CUP,  6,18);  printf)  =========  ==:==========:=====:===::=== 

printf  (CUP,  7,18):  printf)  L  Enter  LI'.'E  caaera  aode .  "  i 

printf  (CUT,  8.18):  printf  (  T  Oisplav  T.iRCA  image  on  monitor  "  i 

'printf(CUP,  9,18);  printf)  '  and  spatial  light  modulator  ") 

printf(CUP,  10 , 18 ); printf (  '  (  .  tga  file  from  disk).  " ) 

printf(CUP.  1 1 ,  1 8  );  printf  :  '  Full  512.v400  .tga  image  is  i 

printf  (CUP.  1 2  ,  18  I  ;  printf  (  "  displayed  on  128x128  SL.'i.  "  ) 

printf  (CUP,  1  3  ,  1 8  )  ;  pr  int  f  (  ■■  0  DISPLAY  only  center  256x255  '  i 

printf(CUP,  1  4  ,  1 8  )  :  pr  intf  (  ■  of  the  laage  on  tne-  128x128  SL.M  '  i 

pri.ntf  (  CUP.  15  .  18  );  printf  ( '■  ') 

printf  (CUP,  16 , 18 1 ; printf (  '  " ) 

printf  (CUP,  I T ,  18 )  ;  printf  (  ■■  Q  Exit  program  (QUIT),  '  i 

printf(CUP,  18  ,  18  );  printf )  ■■  " ) 

printf  (CUP,  19 , 18  );  printf  (  ■■  --press  key  for  desired  function  "  1 

printf (CUP.  20 , 18 ); printf ( "  followed  by  the  ENTER  key.  ") 

printf ( CUP .  23,1); 


keysense ( ) ; 


keys^nse  (  ) 

( 

1 n  t  key; 
int  Found; 

1  r.  t  i ; 

char  inf i ie  f  30 1 : 

iraph Ini t  -  1 i  ; 
ke'-  =  ietchar  (  !  : 

»  w  I  r.  c  n  I IV  e  V  ; 

case  ■ 1 ' : 
case  ’ L ’ : 

Found  =  2  ; 

SetLi veMode (  )  : 
break ; 
case  ' t ' : 
rase  ' T ' ; 

Found  =  1  ; 

printf) "\033[ 40H" ’ ; printf) CLS ) ; printf)  "\033r32H'  ) ; 

pr intf ( CUP . 1 1 . 20 );  printf ( "This  section  will  take  a  .TGA  image  file  and" 
pr  in  t  f  (  CUP ,  1  2 . 20  I  :  pr  int  f  (  display  the  binar'r  version  on  the  Litton  SL.M. 
pr int f ( CUP , 1 5 . 20 1 : pr in t f (  "Please  input  the  TARCA  file  ->  "i; 
scanf)  'Xs\n’',infiiei, 

printf)CUP.20.24  )  (printf  CDISPLAVIMC  I.MACE  D.N  MONITOR  i: 

GetPic) infile. -1. -I, -1 )  : 

printf) CLS ) : 

ccdslmi infile,-l.-l,-li; 
break ; 
case  'd': 
case  'D' : 

Found  =  1  ; 

printf)  "XOSSf  40H"  i  ;  printf)  CLS)  (printf  (  ’  \03  3  [  32H''  )  : 

printf ( CUP , 1 1 . 20  ); printf (  "Only  center  256x256  of  a  .TGA  image  file  will 
printf  (  CUP  ,  12 , 20  );  printf  (  "be  displaved  on  the  128x128  Litton  SL-M."): 


pr intf ( CUP , 1 5 , 20 ) ;printf ( "Pleaae  input  the  TABGA  file  ~>  "l; 
Bcanf ( “XaVa" . infile  I ; 


printf(CUP.20,24l ;printf("DISPLAyiNG  IMAGE  ON  MONITOR" i : 
0etPic( infile,  -1 , -1 ,-l ) : 

printf ( CLS )  : 

ccd«le( infiie,-l,-l,-l); 


«!■(  )  ; 
bre&K ; 
case  ' q ’ : 
case  'O' : 

CraphEndt ) ; 
printf ( “\033 f Oe" ) ; 
printf ( CIS  1  ; 
ex  1 1 1  )  ; 
break ; 
default  : 

Found  =  0 ; 

printf("\033t44; 30a" )  ; 

printf (CUP,  23 . 1) ; printf ( "\033 [ a" ); printf ( "  KEY  NOT  ASSIGNED  i 
for  ii=0;  i<=32000;  i++) 

{ 

/*  pause  to  display  aessage  */ 

! 

printf  ( '■\033  [u"  );  printf  ( '^033  f  K"  I  ;  /»  Erase  line  «/ 
break ; 

I 

if  (Found  ::  0)  keyaense();  /* 

*  Invalid  input.  Get  another  Key. 

*/ 

if  (Found  =s  1)  Benu(  l  ; 
if  (Found  2}  Li'e.Mode(  ) ; 


L  i  l  e.Mode  (  I 
I 

pr intf  (■' \033  (  4  Im"  1  ;  printf  (  CLS  )  ;  /»  Bed  background;  clear  screen  «  / 

pr intf ( "\033 ( 44 ; 32b ■ ) ;  /*  Black  window,  green  letters  */ 


printf (CUP, 
printf ( CUP, 
printf ( CUP , 
pr int  f ( CUP , 
pr intf ( CUP , 
pr intf ( CUP , 
pr intf ( CUP , 
printf (CUP, 
printf (CUP, 
printf ( CUP , 
printf ( CUP, 
pr intf ( CUP , 
printf (CUP, 
printf ( CUP , 
printf ( CLP, 
printf (CUP, 
printf ( CUP, 
printf (CUP, 

Erintf ( CUP, 
ive(  )  : 


3,181  ;  printf  (  " 

4.18) ;  printf ( " 

5.18) ;  printf ( " 

6.18) ;  printf(' 
",18i;  printf('’ 

8.18) ;  printf  ( '■ 

9.18) ;  printf(" 

10.18)  ; printf (  ' 

1 1 . 18 )  ; printf ( " 

12.18)  ;printf ( " 
13,  18) ;printf( " 

14.18)  ;printf (  ’ 

15 . 18 )  ; printf (  ’ 

16 . 18 ) ; printf (  ’ 
ir,18);printf(" 

18.18)  ; printf ( " 

19.18) ;printf(" 

20.18)  ; printf ( " 
23,1); 


LIVE  CAMERA  MODE 
key  function 


C  C.iPTURE  a  video  frame  in  me 

display  aenory. 

S  STORE  image  displayed  on 

monitor  in  a  disk  file. 


M  Return  to  M.ilN  menu. 

Q  Exit  program  (OUIT). 

— press  key  for  desired  function 
followed  by  the  ENTER  key. 
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Live ( ) 
{ 


int  kev; 
int  Found; 
int  i  ; 

char  inf ile  f  50 1 : 

ke V  :  eetchar i i ; 
switch  (key: 

case  ' c  ’  ; 
case  'C: 

Found  =  2 ; 

GrabF  rane (  i  ; 

Graph Init(-ll; 
breuK ; 
case  '  a  '  : 
case  'S': 

Found  =  2: 

printf( "\033(40H" 1 ;printf(CLS» ;printf( ■\033(32H' ) : 
printf ( CUP . 12 , 20  I :printf I "Input  filenaae  ->  " ) ; 
scanf  ( ■'  Xa\n"  ,  inf  ile  )  ; 

print f ( CUP , 1 5 , 20 ) ; pr int f ( "Saving  image  to  disk  as  a  .TGA  file."); 
PutPici infile, -1,-1  ,-l) ; 

pr int f ( CUP , 1 8 , 20 ) ; pr intf ( "Image  saved.  Hit  ENTER  to  continue."); 
getcharf ) ; 
getchar( ) ; 
break : 
case  ' a ' : 
case  ' M ' : 

Found  =  1 ; 
break : 
case  ' q ' : 
case  ’ Q ' : 

GraphEnd ( i ; 
print f :  \03  3  f  Om  '  )  ; 
or int f , CLS  >  : 
exit i  )  ; 
break : 
default 

Found  =  0; 

printfl  ■■\033[44;30m  )  ; 

printf(CUP,  23 , 1 ); print f !" \033 ( s" ); print f i "  KEY  SOT  ASSIGNED"); 
for  (1=0;  i<=32000;  i**i 


/*  pause  to  display 


pr intf ( "\033 [ u" ) ; printf (  XOSO [ K" ) ;  /*  Erase  line  »/ 
break ; 


if  I  Found  ==  0)  Live(); 


if  (Found  =s  II  menul): 
if  (Found  ==  2)  Live^odeii; 


•  Invalid  input.  Get  another  kev. 

•/ 


k  'J^/0  ^ 


.  -• 

*  J*  -  •. 


/ 


Prca:  CCDSL.'t.C  Writtaa  by  Michael  W.  Mayo,  10  Aus  1987 

Cal line  Prograa:  TARGASLM.C 


<  This  function  is  called  by  TARGASLM  to  display  a  TARCA  laase 

*  (.TGAI  onto  the  Litton  spatial  light  aodulator.  The  .TGA  laage 

*  IS  a  312x400  resolution  image;  each  pixel  having  an  3  bit  grav 

*  scale  value  between  0  and  233.  This  function  takes  the  312x400 

<  laage  and  ..'oapresses  it  to  a  128x123  binary  image  (a  pixel  on 

*  on  the  SLM  is  either  on  or  off  -  no  gray  levels). 

t 

t  CCDSLM.C  IS  a  aodified  version  of  the  TARGA  library  function 
t  GETPIC.C.  Instead  of  opening  the  .TGA  disk  file  and  dumping 
«  it  to  the  nigh  resolution  monitor,  the  file  is  intercepted, 

*  compressed,  binarized,  and  sent  to  the  SL.M. 

t 
t  / 


^include  <stdio.h> 

•include  tardev.h" 

Sdefine  MAXBUFSI2E  32767  /»  must  be  at  least  2048  */ 

•define  0N_MASK  0x80  /»  Allows  write  to  SU*  by  putting  SPA7  high 

•define  OFF_MASK  0x7f 

extern  char  *aalloc(  )  ; 

extern  struct  TARStruct  *targa; 
static  struct  hdStruct  header; 
extern  char  *strehr(): 

extern  unsigned  char  pokeb( 1 ; 

unsigned  legmentsOxbOOO ;  /«  Offset  segment  of  poke  address  »/ 

unsigned  simrow.  slmcol;  /*  Ro-  and  col  of  SL.'I  .0-127; 

char  byte;  /*  Bvte  data  sent  to  poke  address  «' 

int  ccdsimi picture .XCorner,  VCorner.  screen) 

char  picturef 1 ; 

int  XCorner,  VCorner; 

int  screen; 

/ 

char  ♦buffer; 
char  *CoiorMap; 
char  »dataPntr: 
int  pict ; 

unsigned  buffSize; 

ch.ir  inname  t  30]; 

int  points; 

int  xOffaet,  ''Offset; 

int  RowaPerRead  ,  ByteaPerRead ,  ByteslnPi.xel : 

ir.t  BytesInRow; 

int  row,  1,  col: 

int  mapBytea; 

int  mapLcngth; 

int  Convert:  /*  flag  to  indicate  type  of  convection  required  •/ 

int  Compressed;  /*  flag  to  indicate  if  Compressed  image  */ 

int  InputBytes;  /»  number  o  bytes  per  input  pixel  */ 

int  PixPerRead;  /*  number  of  pixels  in  one  read  */ 

int  intensity;  /*  0-255  intensity  value  of  a  pixel  */ 

int  threshold;  /*  Threshold  value  for  binary  write  to  SLM  */ 


i:-:" 


/  • 

t  GET  T\BCA  FILE 

t 

<  Coapil*  the  naae  first 
»/ 


strcpT ( innaae . picture ) ; 

if  (  strchr  (  innaae  ,  '  .  '  )  =  =  ML'LL '•  street  I  innaae  .  ~  .  tsa\0"  )  ; 


/» 

*  .  7pen  the  file  header  and  see  if  a  valid  tarca  iaage 
»/ 


I  f  (  CetPicHeader (  innaae  .lipict , iiheader ,  0  )  ==  -ll  I 

printf ( "XnVr  IMAGE  Xs  NOT  FOUND\n\r  ' , picture ) : 

return ( -  1  )  :  ' 

) 

BytesInPixel  =  tarja- > By tesPerPixel ;  /*  bytes/pixei  in  display  «/ 

if  (BytesInPixel  ==  3  t  BytesInPixel  =4: 

InputBytess  ( header . dataBita+T ) /8 ;  /*  byte/pixel  on  the  disk  */ 


/* 

*  DEFINE  IMAGE  ORIGIN 

« 

«  set  the  corner  unless  specif iec  otherwise 
*/ 

i.  XCorner  ==  -1  ) 

ixOffset  s  header . XOf f set :  yOffset  =  header .YOffset :  1 

else 

I  xOffset  =  XCorner;  yOffset  s  YCorner:  I 


Converts  0;  /»  default  fla<  no  conversion  necessary  »/ 

/*  .f  color  mapped  covert  the  map  as  needed  */ 


if  1  BytesInPixel  '=  InputBytes  )  Convert  -  1; 


/* 

*  See  if  compressed  file  or  not 

*/ 

if  I  header . datatype  >  8  )  Compressed  =  1; 
else  Compressed  =  0; 

/* 

»  SET  PARAMETERS  TO  ALLOW  MY  DISK  READ 

»/ 


/* 

« 

t 

« 


points  =  header. x: 

Compute  the  maxiaua  number  of  rows  that  we 
can  read  in  at  a  shot  and  the  number  of  bytes 
to  read  from  the  disk  (assume  uncompressed). 


if  (  BytesInPixel 
else 


InputBytes  ) 


( Bytes InRow  = 
{ Bytes InRow  = 


BytesInPi.xel*points ;  1 
InputByte8»points ;  ) 


if 


(Compressed  ==  0  )  ( 

buffs  i-ze  =  .MA.XBUFSIZE. 

while  (  (  (  ( int ) ( buf fer=maiioc ( buf fS ize  I  I 

( buffs ize> sBytes InRow I )  { 

buffSize  =  (  buf  f  Size'*’ 1  )  >  >  I : 

) 

if  (  buffSize  <  BytesInRow  )  { 

close ( pict ) ; 
return( -2  I  ; 

I 


.NIL'LL)  lA 
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)  ( 

buffSize  s  B7t.eBlnKow: 
buffer  =  Belloc(  BjteelnSow); 
if  (  ExpStarK plot , InputBytee )  <  0  )  { 
free (buffer ) ; 
closet  pict ) ; 
return( -2 ) ; 


RowsPerRead  =  buffSise/BytesInRow; 

ByteaPerRead:  RowsPerRead  *  InputBytea  t  points; 
PixPerRead  =  R^wiPerRead  (points ; 

/* 

*  Erase  the  SLM  display 

»/ 

erase i i ; 

/» 

*  Targa  iaage  sent  to  spatial  light  modulator. 

»  / 

printf { "\033(  15:22H" )  ; 

pr intf (" Input  isage  intensity  threshold  value:  '); 
scanfl  "  Xd  \n".  ithreshold):  printf  ( ■■\033 1  2  J"  )  ; 
printf ( "\033[ 10;26H’ ) ; 

printf ' "Transi erring  iaage  to  the  SLM.'); 

slarow  =0:  /»  First  row  on  SLM  */ 

for  I  rowsO;  row<  header,  y;  row  •♦•=  RowsPerRead  ) 


*  Read  the  date  froa  the  TARGA  file 

*/ 

if  (Compressed  =s  0)  read( pict , buffer , BytesPerRead ) ; 
else  ExpRead( buffer. PixPerRead) ; 

i f  ( Convert  =s  1 ) 

LnPacitBuf  (  buffer  .ColorMap.Pi.vPerRead.BytesInPisei )  : 
dataPntr  =  buffer; 

for  (1=0;  ( i tRowsPerRead )4t( { itrow i <=( header . >-3 )) :  !♦♦) 


*  Convert  TARGA  iaage  froa  0-255  gray  levels  to  binaray  and 

*  place  binary  iaage  on  SLM. 

*/ 

if  ( ( rowfi ) <=7 )  dataPntry=512 :  /» 

*  Skip  1st  8  rows  but 

*  update  pointer. 

»/ 

else 

( 

if  ( ( ( row+i+l )X3 )=s0)  /t  Take  every  3rd  row.  »/ 

( 

s1bco1=0 ; 

for  (col=0;  col<=511;  col+y ) 


if  ( ( (col)X41==0» 

( 


/»  Every  4th  coluan  »/ 


*  Cray  levels  are  at  -128  to  127.  Convert  so  they  are  0  to  255 

«/ 

intens ity= (dataPntr ; 

if  (  intensitytO  )  intensity'^  =  256 ; 

if  ( inten8ity>threshold) 


byte  =  slarow  ;  ON_'1ASK;  Assures  SPA7  is  high  */ 

pokeb  (  slacol .  segaent .  by  tc  )  ;  /»  Write  pi.yel  to  SLM  */ 


/<  Nazi  col  on  SLM  </ 


» 


•  lacol^'f: 

) 

dntnPntr^-*- :  /•  Update  pointer  */ 

) 

alnrow'^'f;  /»  Next  row  on  SLM  »/ 

) 

else 

dataPntr  ♦=  512:  /»  Skip  row  but  update  pointer  */ 

I 

,  I 

) 

pokeb( 2049 , aegnent , 255 ) ;  /*  Coil  pulae  to  write  inage  on  SLM  */ 

printf ( "\033[ 19;27H" )  ; 

printf ( " laage  dump  to  SLM  coaplete."); 
pnntf  (  'AOSOI  24  :23H"  )  ; 

printf ( "Preaa  ENTER  to  return  to  Main  Menu"); 
getchar( );  /»  Waiting  for  ENTER  key  */ 

getchar( ) ; 

/* 

*  Write  to  SLM  coaplete. 

</ 

if  (  Conpreaeed  ==  1  )  ExpStop( ) : 
f  ree ( buf  f  er )  ; 
ciose'pict) ; 
return ( 0  I  ; 


/* 

<  This  function  ia  used  to  clear  (erase)  the  spatial 
*  light  nodulator  display. 

»/ 


erase (  ) 

I 

int  row,  col; 


) 


printf( "\033(5;30H"  )  ; 

printf ( "Erasing  SLM  display."); 

for  (row=0;  row<  =  l27;  row'*-'f  ( 

( 

for  (col=0;  col<=127;  col^+l 

( 

byte  =  row; 

pokeb( col , segment , byte  ) ; 

1 

» 

pokeb( 2049 , segment , 255  )  ; 

printf("\033[8;28H")  ; 

printf ( "Erase  operation  complete. 


/*  Places  cursor  at  »/ 
/»  row  5  and  col  30  »/ 


/» 

Nucleate 

each 

row 

*/ 

/* 

and  column  on 

display 

»/ 

/* 

Saturate 

( coil 

pulse  ) 

to 

complete 

erase 

»/ 

75 


/•  Pm:  SLM.C  Written  by  Michael  W.  Mayo.  5  Sep  1987 

t  Called  bT  TARGASLH.C 

t 

«  This  prograa  is  very  siailar  to  CCDSLM.C.  A  .TCA  iaage 

*  file  is  retreived  froa  disk,  compressed,  binarized,  and  sent 

*  to  the  spatial  light  aodulator.  what  makes  this  function 

*  different  is  that  only  the  middle  256x256  of  the  512x400  image 

*  is  used  (in  other  words,  the  outside  of  the  image  is  cut  off). 

<  The  remaining  256x256  image  is  then  sampled  (every  other  pixel) 

<  and  displayed  as  a  128x128  image  on  the  spatial  light  aodulator. 

*  This  program  is  needed  to  increase  the  image  resolution  on 
t  the  spatial  light  aodulator  and  decrease  the  sampling  error. 

*  SLM.C  was  used  mainly  to  display  Fourier  transforms  of  images  on 
»  the  SL-M. 


sinclude  <stdio.h> 
•include  "tardev.h" 


idefine  HAXBLFSIZE  32767 
•define  O.V_MAbK  0x80 
•define  OFF  MASK  0x7 f 


/*  must  be  at  least  2048  */ 

/*  Allows  write  to  SLM  by  putting  SPAT  high  */ 


extern  cnar  *malloc(); 

extern  struct  TARStruct  *targa: 
static  struct  hdStruct  header; 
extern  char  tatrchr( ) ; 


extern  unsigned  char  pokeb( ) ; 
unsigned  segaentsOxbOOO ; 
unsigned  slarow,  slacol ; 

Char  byte; 


/»  Offset  segment  of  poke  address  »/ 
/»  Row  and  col  of  SLM  (0-127)  »/ 
/»  3yte  data  sent  to  poke  address  */ 


int  slm( picture . XCorner ,  YCorner.  screen) 

char  picture ( 1 : 

int  XCerner,  YCorner; 

int  screen; 

( 

char  ^buffer; 
cha,'-  *ColorMap; 
char  »dataPntr; 
int  pict: 

unsigned  buf f Size ; 


char  innane[ 50 ) ; 

int  points; 

int  xOffset,  yOffset; 

int  RowsPerRead,  ByteaPerRead ,  Bytes InPixel ; 

int  B*  eslnRow; 

int  r<  i  .  col ; 

int  ar.:  'tes ; 

int  ma,.wengtb; 

int  Convert;  /*  flag  to  indicate  type  < 

int  Coa:iressed;  /*  flag  to  indicate  if  Coi 


int  InputBytes; 
int  PixPerRead; 
int  intensity; 
int  threshold; 


flag  to  indicate  type  of  convention  required  */ 
flag  to  indicate  if  Compressed  image  */ 
number  o  bytes  per  input  pixel  */ 
number  of  pixels  in  one  read  */ 

0-255  intensity  value  of  a  pixel  */ 

Threshold  value  for  binary  write  to  SLM  »/ 


GBT  TABCA  FILB 


« 

t 

* 

«/ 


Coapile  tha  naa«  first 


strcpjr^  innaae , picture  i  ; 

if  (  Btrchr  (  innaae ,  '.'>=s!4ULL)  strcat  (  innaae .  "  .  tcaVO"  I  ; 


< 

*  open  the  file  header  and  see  if  a  valid  tar(a  iaase 
»/ 


if (GetPicRpader(innaae,ftpict,&header.O)  =:  -II  ( 

printf  (  "XnVr  IMAGE  Xa  NOT  FOUNDynSr** , picture  I ; 
return(-ll; 

) 

BytealnPixel  =  targa- >BTteaPerPixel ;  /»  bytea/pixel  in  display  */ 

if  (BytealnPixel  ==  3  I  BytealnPixel  =4; 

InputBytea:  (  header  .dataBitS'^7  ) /8 ;  /*  byte/pixel  on  the  disk  */ 


<  DEFINE  IMAGE  ORIGIN 

« 

*  set  the  corner  unless  specified  otherwise 
» / 

if  (  XCorner  ==  -1 

(xOffaet  =  header .XOff set;  yOffaet  =  header . YOffaet ;  ) 
else 

(  xOffset  s  XCorner;  yOffaet  =  YCorner;  ) 

Converts  0;  /*  default  flag  no  conversion  necessary  </ 

/*  if  color  aapped  covert  the  asp  as  needed  */ 

if  (  BytealnPixel  InputBytes)  Convert  =  1: 

< 

»  See  if  coapresaed  file  or  not 

»/ 

if  (  header .dataType  >  8  )  Coapresaed  s  i; 
else  Coapresaed  s  0; 

t 

*  SET  PARAMETERS  TO  vLLOW  MY  DISK  READ 

*/ 

points  s  header. x; 

I 

*  Coapute  the  aaxiaun  nuaber  of  rows  that  we 

t  can  read  in  at  a  shot  and  the  nuaber  of  bytes 

*  to  read  froa  the  disk  (assuae  uncoapresaed ) . 

«/ 

if  (BytealnPixel  >  InputBytes)  (BytealnRow  =  BytealnPixel *points ; ) 
else  (BytealnRow  =  InputByteatpointa ; ) 

if  (Coapresaed  ==  0  )  { 
buffSize  =  MAXBUFSJZE; 

while  (  (  (  ( int )( buffer=aalloc  (buffSize  I  I  )  -=  .'.'ULD  && 
(buffSize>=BytesInRow) )  ( 

buffSize  =  ( buf fSizee 1 ) > > 1 ; 

) 

if  (  buffSize  <  BytealnRow  )  ( 

close) pict) ; 
return) -2  I : 


f 

Ad 


I 


♦s 

I 

•s 


'I 


'a 

.♦S 


1  «l»e  { 

buffSize  s  BrtesInRow; 
buffer  3  m11oc(  BjrteelnRow) ; 
if  (  ExpStart ( pict , InputBytes )  <  0  )  ( 

free (buffer ) ; 
close ( pict I : 
return ( -2  )  ; 


RowePerRead  =  buf  fSize/ByteslnRow ; 

BytesPerRead:  RowaPerRead  <  InputBytes  *  points; 
PixPerRead  :  RowaPerRead  <points; 

/» 

<  Erase  the  SLM  display 
»  / 

erase i ' ; 

/» 

*  Tarsa  image  sent  to  spatial  light  modulator. 

»/ 

printf (  AOaai 15:22H" ) ; 

printt(  Input  image  intensity  threshold  value:  "); 
scanf (  "  \d  \n" ,  ithreshold);  printf ( "NOSS ( 2 J" ) ; 
printf ( "NOaaC 10;26H" ) ; 

printf I "Transferring  image  to  the  SLM.">; 

slmrow  =  0;  /»  First  row  on  SLM  »/ 

for  (  row=0:  row<  header. y;  row  »=  RowaPerRead  ) 


*  Read  the  date  from  the  TARCA  file 

»/ 

if  (Compressed  ==  0)  read( pict , buffer , BytesPerRead ; : 
else  ExpRead ( buffer (PixPerRead) ; 

1  f  ! Convert  ==  1 ) 

CnPackBuf ( buffer , Co 1 orMap . PixPerRead .BytesInPixel: ; 
aataPntr  =  buffer: 

for  (1  =  0;  I  1 <RowsPerRead i a& I  I itrow I < = t header . v-i 2 ) I :  i*+) 

( 

' « 

*  Conver*^  TARGA  image  from  0-2S5  to  binaray  and 
<  place  binary  image  on  SLM 

*/ 

if  ( ( rowt i ) <  =  7 1 )  dataPntrtsSlE;  /» 

«  Skip  Ist  72  rows  but 
*  .update  pointer. 

»/ 

else 

( 

if  ( ( ( rowt i )X2 )=s0)  /*  Take  every  other  row.  «/ 

{ 

slmcolsO ; 

dataPntrt=128;  /•  Skip  first  128  cols  »/ 

for  (col=0;  col<=255;  col++) 

( 

if  ( ( (coI)X2)==0)  /*  Every  other  column  »/ 

( 

intensity3*dataPntr ; 

if  (  intensitytO)  intensity't-3256 ; 

if  ( intensity>threshold) 

( 

byte  =  slmrow  ;  ON_MASK: 

pokeb( slmcol . segment . byte ) ;  /*  Write  pixel  to  SLM  */ 


/*  Update  pointer  */ 

/*  Skip  last  12S  columns  */ 


bIbcoI^^'*'  ; 

) 

dataPntr+^ ; 

) 

dataPntr*=128 ; 
slarow«« ; 

else 

dataPntr  4-=  312:  /*  Skip  row  but  update  pointer  */ 

) 

» 

t 

pokeb( 2049 , aecaent , 255 ) :  /•  Coil  pulse  to  write  image  on  SLM  */ 

printf{"\033(19:27H"); 

printf (” Image  dump  to  SLM  complete."); 
printf ( "\033(24;23H" ) ; 

printf ( "Press  ENTER  to  return  to  Main  Menu"); 
getchar)  );  /*  Waiting  for  ENTER  kejr  */ 

getchar( ) ; 

t 

*  Write  to  SL.M  complete. 

*/ 

if  t  Compressed  ==  1  )  ExpStopt ) ; 
f  ree ( buf  f er i ; 
close ( pic  1 1 ; 
return ( 0 ! ; 


erase()  /<  Erases  SLM  display  */ 

1 

int  row,  col; 

printf ("\033[5:30H") ; 

printf  (  "Erasing  SL.M  display."); 

for  irowsO;  rowCiZT;  row+t! 

r 

:'or  (ooi:0:  coi'.  =  12T;  col't’+l 
) 

byte  =  row; 

pokeb( col , segment ,byte ) : 

) 

I 

pokeb( 2049 , segment , 255  )  ; 
printf) "\033[8;28H" )  ; 

printf ( "Erase  operation  complete."!; 
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Fabrication  of  the  GGH 


This  Appendix  describes  the  proceedure  of  producing  and 
printing  the  grating  pattern  which  will  perform  the  desired 
optical  coordinate  transformation  (CT) .  Also  included  is 
the  method  used  to  reduce  and  place  the  grating  pattern  onto 
a  high  resolution  film  plate. 

Computer  Generated  Intef erogram 

The  computer  generated  inteferogram  is  made  up  of  a 
binary  fringe  pattern  printed  on  a  laser  printer  and  then 
photoreduced .  Due  to  the  binary  nature  of  the  intefero¬ 
gram,  an  implicit  carrier  frequency  is  inherent  in  the 
hologram.  As  a  result,  many  diffracted  waves  will  result 
when  the  CGH  is  illuminated.  A  sufficiently  high  carrier 
frequency  must  be  chosen  to  seperate  the  first-order  dif¬ 
fracted  wave  (where  the  CT  is  located)  from  the  higher-order 
waves.  To  avoid  aliasing  (overlapping  of  diffracted 
orders) ,  the  following  equation  must  be  satisfied  for  the 
carrier  frequency  a  (6:  939) 

°  ]  ""[  (  I  %.x]'  *  (  5  w)' 

where:  and  y  are  input  dimensions, 

max  max 

^  is  the  wavelength  of  light  used, 
f^  is  the  focal  length  of  CT  lens. 


The  CGH  constructed  for  this  thesis  effort  had  the  following 


parameters:  x  =  y  =  10mm  (this  is  the  dimensions  of 
max  max 

the  spatial  light  modulator);  i  =  0.6328^m;  and  f.  =  300mm. 

Li 

Using  these  parameters  in  Equation  B.l  and  solving  results 
in  an  a  >  30.9  lines/mm.  An  a  of  35  lines/mm  was  chosen  to 
meet  the  above  criteria. 

To  form  the  computer  generated  hologram  (interfer- 
ogram) ,  the  binary  fringes  of  the  transmittance  function 
were  printed  on  a  laser  printer  and  then  optically  reduced. 
For  the  lnr-0  CT  CGH  the  transmittance  function  obtained  by 
interfering  the  required  phase  function  given  by 


0  (x,  y ) 


=  r  —  1* 


,  ,  2  ^  2,1/2  ^  ,  /  V 

ln(x  +  y  )  -  ytan(y/x) 


-X 


(8) 


with  a  carrier,  a,  is  given  by 

t(x,y)  =  O.sj^l  +  cos[2rtax  -  0(x,y)]j  (B.2) 

The  transmittance  function  has  fringes  (maxima)  at  positions 
(x,y)  satisfying 


2^ax  -  0(x,y)  =  2nnj^  (B.3) 

for  various  values  of  where  each  value  of  corresponds 

to  a  different  fringe. 

The  various  (x,y)  points  that  satisfy  (B.3)  for  each 
value  of  nj^  were  connected  to  form  the  fringe  pattern. 

The  Metalib  program  (for  contours)  on  the  AFIT  computer 
systems  was  used  to  solve  for  the  x,y  points  and  the 
resulting  fringe  pattern  was  printed  on  the  Imagen  laser 
printer.  The  CGH  plot  is  shown  in  Figure  B.l. 


After  the  plot  was  obtained  from  the  Imagen  Laser 


i  Printer,  a  negative  of  the  plot  was  produced  at  the  Base 

Photo  Lab.  The  negative  was  taken  to  the  Cooperative 
Electronics  Materials  Processes  Lab  (room  1065)  in  building 
125  to  be  reduced  onto  a  high  resolution  film  plate. 

The  Dekacon  optical  system  was  used  to  reduce  the  CGH 
and  expose  the  film  plate.  A  3-inch  Wray  lens  was  placed 
into  the  system,  as  shown  in  Figure  B.2,  for  a  20x 
reduction.  A  20x  reduction  was  obtained  by  setting  the 
front  box  at  6.883  and  the  rear  box  at  27.000.  The 
resulting  image  after  these  settings  was  lOxlOmm  with  a 
carrier  frequency  (a)  of  35  lines/nun. 

The  negative  of  the  CGH  plot  was  placed  on  the  object 
screen  and  back  illuminated  with  green  light.  Red  Rubylith’^ 
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Figure  B.2.  Optical  reduction  of  CGH. 
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material  was  placed  around  the  plot  to  block  all  remaining 
light;  therefore,  the  only  light  from  the  object  screen  was 
that  which  passed  through  the  CGH  plot. 

Scratches  were  placed  on  a  clear  plate  and  placed  in 
the  image  plane  of  the  system.  A  microscope,  placed  behind 
the  image  plane,  was  used  to  focus  on  the  scratches  residing 
in  the  image  plane  while  the  front  box  was  moved  until  the 
CGH  plot  came  into  focus  with  the  scratches.  When  the  plot 
and  scratches  were  both  in  focus,  the  system  was  focused  for 
the  desired  reduction.  Note:  it  is  very  important  that  the 
scratched  surface  of  the  glass  faces  the  object  screen  since 
the  emulsion  on  the  film  plate  will  also  face  that 
direction , 

Each  Kodak  2x2  in.  High  Resolution  Plate  was  exposed 
for  60sec  by  removing  the  shutter  from  the  plate  holder. 

The  exposed  plates  were  then  developed  for  40-60sec, 
depending  on  the  type,  age,  and  temperature  of  the 
developer.  Best  results  were  obtained  by  pulling  the  plate 
out  of  the  developer  every  few  seconds  (after  30sec)  and 
viewing  it  under  red  light  until  the  grating  pattern  seemed 
"dark  enough".  The  developed  plate  was  then  placed  in  the 
stop  bath  for  30sec  followed  by  one  minute  in  the  fixer.  A 
five  minute  rinsing  in  deionized  water  completed  the 
development  process.  The  developed  plate  was  then  viewed 
under  a  microscope  to  assure  proper  results  (lines  should  be 
very  dark  and  spacing  between  lines  should  be  very  clear) . 


Bibliography 


L.  "AT&T  Truevision  Advanced  Raster  Graphics  Adapter; 

Targa  8;  Software  Tools  Notebook,"  operation  manual  for 
the  frame  grabber  system:  AT&T  Electronic  Photography 
and  Imaging  Center,  Release  3.0  (November  1986). 

2.  Casasent,  David.  "Coherent  optical  pattern 
recognition:  a  review,"  Optical  Engineering,  24: 

26-32  (January  1985) . 

3.  Casasent,  David  and  A.  Furman.  "Sources  of  Correlation 
Degradation,"  Applied  Optics,  16:  1652  (1977). 

I.  Casasent,  David  and  Hironob  Okuyama .  "A  High- 
Dimensionality  Pattern  Recognition  Feature  Space," 

SPIE  Intelligent  Robots  and  Computer  Vision,  579: 
245-256  (1985) . 

5.  Casasent,  David  and  D.  Psaltis.  "Deformation 

Invariant,  Space-Variant  Optical  Pattern  Recognition," 
Progress  in  Optics,  16;  289  (1978). 

5.  Casasent,  David  and  others.  "Real-time  deformation 

invariant  optical  pattern  recognition  using  coordinate 
transformation,"  Applied  Optics,  26;  938-942  (March 
1987) . 

7.  Cederquist,  J.  and  A.  Tai.  "Computer  Generated 

Holograms  for  Geometric  Transformations,"  Applied 
Optics,  23;  (18)  3099-3104  (September  1984). 

8.  Flannary,  David  and  others.  "Real-time  coherent 

correlator  using  binary  magneto-optic  spatial  light 
modulators  at  input  and  Fourier  planes,"  Applied 
Optics,  25:  (4)  466  (February  1986). 

9.  "Functional  Description  for  the  Low  Current  Driver  Card 
Set  for  the  128x128  Display  System  for  the  L135 
Program,"  operation  manual  for  the  Litton  Data  Systems 
128x128  LIGHT  MOD. 

10.  Goodman,  J.W.  Introduction  to  Fourier  Optics.  New 
York:  McGraw-Hill  Book  Company,  1986. 

II.  Gregory,  Don  A.  "Real-time  pattern  recognition  using  a 

modified  liquid  crystal  television  in  a  coherent 
optical  correlator,"  Applied  Optics,  25:  (4)  467-469 

(February  1986) . 


85 


12.  Hughes,  Kenneth  D.  Suitability  and  Applications  of 

Liquid  Crystal  Televisions  In  Optical  Pre-processors. 
MS  thesis,  AFIT/GE/ENG/86D-6 .  School  of  Engineering, 
Air  Force  Institute  of  Technology  (AU) ,  Wright 
Patterson  AFB  OH,  December  1986. 


i 


li 


13.  Jensen,  A.S.  and  others.  "Transformation  of  image 
positions,  rotations,  and  sizes  into  shift  parameters," 
Applied  Optics,  26:  (9)  1775-1781  (1  May  1987). 

14.  Kobel,  W.G.  and  T.  Martin.  A  Distortion-Invariant 
Pattern  Recognition  AlgorithiiK  MS  thesis, 
AFIT/GE/ENG/86D-20 .  School  of  Engineering,  Air  Force 
Institute  of  Technology  (AU) ,  Wright-Patterson  AFB  OH, 
December  1986. 

15.  Neidig,  David  L.  Optical  Pattern  Recognition  Using 
Synthetic  Discriminant  Functions,  MS  thesis,  GE/EE/86D. 
School  of  Engineering,  Air  Force  Institute  of  Tech¬ 
nology  (AU) ,  Wright-Patterson  AFB  OH,  December  1986. 

16.  Psaltis,  Demetri  and  others.  "Optical  image 
correlation  with  a  binary  spatial  light  modulator," 
Optical  Engineering,  23;  698-704  (December  1986)  . 

17.  Ross,  W.E,  and  others.  "Fundamental  characteristics  of 
the  Litton  iron  garnet  magneto-optic  spatial  light 
modulator,"  SPIE  Advances  in  Optical  Information 
Processing,  388;  55-64  (1983) . 

18.  "Semetex  Corporation  SIGHT  MOD  Development  System 
Operations  Manual,"  for  the  IBM-PC/SIGHT  MOD  interface 
board . 


o'hjs. 


VITA 


Michael  H.  Mayo  was  born  on  28  August  1964  in 
Wiesbaden,  Germany.  He  graduated  from  O' Fallon  Township 
High  School  in  O' Fallon  Illinois  in  June  1982.  Having 
earned  an  Air  Force  ROTC  scholarship,  he  attended  the 
University  of  Central  Florida  in  Orlando,  Florida, 
graduating  in  May  1986,  with  honor.  He  received  the  degree 
of  Bachelor  of  Science  in  Engineering,  electrical 
engineering  option  with  a  specialty  in  electro-optics.  Upon 
graduation,  he  received  a  commission  in  the  USAF  through  the 
ROTC  program.  After  commissioning,  he  entered  the 
Electro-optics  Masters  Program  in  the  School  of  Engineering, 
Air  Force  Institute  of  Technology,  in  May  1986. 


Permanent  address:  205  Dartmouth  Dr. 


O'Fallon,  II  62269 


UNCLASSIFIED 


IFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


1«.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


SECURITY  CLASSIFICATION  AUTHORITY 


1b.  RESTRICTIVE  MARKINGS 


gjj 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

AF IT /GEO  /ENG  /8  7D- 1 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

School  of  Engineering 


6b.  OFFICE  SYMBOL  I  7a.  NAME  OF  MONITORING  ORGANIZATION 
(If  applicable)  I 


AFIT/ENG 


6c  ADDRESS  {City.  State,  and  ZIP  Code) 

Air  Force  Institute  o£  Technology 
Wright-Patterson  AFB,  Ohio  45433 


7b.  ADDRESS  (C/ty,  State,  and  ZIP  Code) 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


Bb.  OFFICE  SYMBOL  I  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable)  I 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PROJEa 
ELEMENT  NO.  I  NO. 


8c  ADDRESS  (City.  State,  and  ZIP  Code) 


1 1 .  TITLE  (IfKiude  Security  Clasiification) 

CCMPUTER  GENERATED  HOLOGRAM  AND  MAGNETO-OPTIC  SPATIAL  LIGHT  MODULATOR 


12.  PERSONAL  AUTHOR(S) 

_ Michael  W.  Mayo,  B.S. ,  Second  Lieutenant,  USAF 


.  TYPE  OF  REPORT  13b.  TIME  COVERED  1 14.  DATE  OF  REPORT  (Year,  Monf/f,  Day^  15.  PAGE  COUNT 

MS  Thesis  FROM _ to _  I  1987  December  96 


COSATI  CODES 


GROUP 


05 


SUB-GROUP 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

tot  leal  Pattern  Recognition,  Computer  Generated 
Holcjgram,  Spatial  Light  Modulator 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Title;  CCMPUTER  GENERATED  HOLCX3RAM  AND  MAGNETO-CPTIC  SPATIAL 
LIGHT  MODULATOR  FOR  Cf’TICAL  PATTERN  RECOGNITION 


Thesis  Chairman; 


Dr.  Matthew  Kabrisky 

Professor  of  Elec±rlcal  Engineering 


. 

Ml.  m  .nd  Prt.l,„iono|  D«T«|«aMa 

”  cl  IcchacIcffT  [ftlar 

Afa  OH 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  IBI.  ABSTRACT  SICURITY  CLASSIFICATION 

<^^UNeLASSIFIID/UNLlMITSO  □SAMIASRPT.  □  OTIC  USERS  I  UNCLASSIFIED  _ 


name  of  RESPONSIBLE  INDIVIDUAL  |22b.  TELEPHONE  (Include  Area  Code)  1 22e.  OFFICE  SYMBOL 

Dr.  Matthew  Kabrisky 


00  Form  1473.  JUN  86 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAG 


UNCLASSIFIED 


UNCLASSpiED 

Security  Classification  of  This  Page 


This  thesis  investigates  the  integration  of  an  optical 
system  for  real-time  position,  scale,  and  rotation  invariant 
pattern  recognition.  Specifically  a  Litton  Magneto-Optic 
Spatial  Light  Modulator  is  interfaced  to  a  Zenith  248 
microconputer  and  AT&T  frame  grabber.  A  user  interface 
written  in  C  allowed  arbitrary  patterns  to  be  written  to  the 
SLM  by  a  user.  Patterns  can  be  generated  in  computer  code, 
transferred  from  a  CCD  camera  (via  the  frame  grabber) ,  taken 
from  an  image  file  on  computer  disk,  or  taken  fran  a  VAX 
image  file.  The  ability  to  locally  generate  a  computer 
generated  hologram  using  the  computer  generated  interfero- 
grem  method  on  the  Sun  workstations  and  Imagin  laser  printer 
was  developed.  The  CGH  developed  was  used  in  an  optical 
coordinate  transformation  to  create  a  position,  scale,  and 
rotation  Invariant  feature  space.  CCD  cameras  were  used 
throughout  the  optical  processor  for  display /analysis .  The 
entire  system  was  developed  and  tested  using  various 
templates  and  real  scenes. 


UNCLASSIFIED 


